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Abstract 


The objective of this thesis is to deal with linear integro-differential Fredholm equations especially as 
these kinds of equations play an important role in modeling of different problems in various fields such 
as physics and biology. Therefore, we are going to study them in the analytically and numerically sense. 

We study these equations in Banach space C'/a, b] with two cases of continuous and weakly singular 

kernels. In the continuous case we construct three methods based on the Nystrém, Collocation and 

Kantorovich methods in order to find the best approximation of our solution. In the weakly singular 

case, we construct two methods which are b-spline collocation and product integration in an essential 

reason which is a good precision and acceleration in the calculations.. 

We analysis our equation analytically in Sobolev spaces W1?{a, b], p € [1,-+00|. We give a sufficient 
condition that shows the existence and uniqueness of the solution in these mathematical spaces. We 
have constructed in both spaces: W1![a,b] and H*{a,b] two projection methods based on Galerkin and 
Kantorovich. 


Keywords: Fredholm integral equation, Integro-differential equation, System of integral equations, 
Projection methods, Iterative Methods, Nyst6ém Method 


Mathematics Subject Classification: 45B05, 47G20, 45F05, 65R10, 65F10,64R20 


Résumé 


L’objectif de cette these est de traiter les équations linéaires intégro-différentielles de Fredholm 
analytiquement et numériquement, surtout que ce genre d’équations, jouent un réle important pour 
modeéliser différents problémes dans plusieurs domaines, en particulier la physique et la biologie. 
Nous allons étudier ces équations dans l’espace de Banach C'Ja, b], dans les cas respectivement des 
noyaux continues et faiblement singulieres. 

Dans le but de chercher la meilleur approximation, nous allons construire dans le cas continue, trois 
méthodes numériques, basées sur les méthodes Nystrém, Collocation et Kantorovich. 

Dans le cas faiblement singuliére, nous adaptons le principe de comparaison en appliquant deux 
méthodes, la méthode des intégrations de produits et la méthodes b-spline collocation, afin déterminer 
la bonne précision concernant l’erreur et le temps de calcul. 

Nous traitons notre équation analytiquement dans les espaces de Sobolev W1?[a, b], p € [1, +00] et 
nous proposons une condition suffisante, pour montrer l’existence et l’unicité de la solution dans ce 
types d’espace. Enfin, nous adaptions notre techniques aux espaces : W![a,b] et H+[a, b] en associant 
deux méthodes de projection a savoir Galerkin et Kantorovich. 


Mots clé : Equation intégrale de Fredholm, Equation Integro-differential, Systeme des équations 
integrales, Méthodes de projection, Méthodes Itératives, Méthode de Nyst6m 


Mathematics Subject Classification :45B05, 47G20, 45F05, 65R10, 65F10,64R20 


Qaeddtall 


9 BSI! dvaai ly) Ul 9 Aes! AL oa! AotW! elgrty,3 GVoleo ao Saleill ss dag bY! ode (ys Cargll 
Le gloul 9 sbjudll earl yo drolall Meal! alice 9 Slivall yo tell doded (3 Yoltsinl IMS yo LS 
Goda)l Guodly des JSiy inl ty agai Gyo 


ORR E (51939 Qu poiane Quiles Quilgd yurdlicne WIL> ao C1[a, b] TLL clad (8 SVrlroI! ode Guy 
OF (dg 9395 9 Quam! dd b a9 ids | @ Gb SW Ne BSG yo And) gle 48 Brine Ved 
cde aleteVb gins lomgiyadinte 20 pete gutlieas qulledl dle 8 San gla do LB leed Jal 
bluse 3 Acyl 9 BUI LG Go Gad) Uo ail le geod Gag gitall ey b-spline Jig 


Fol Ibm 3 oll Avlings 9929 Wad dudsd ely UME yo WS 9 Ut pg DleLad (3 Wed Lioles Guy 


9 Wt [a,b]: SlelLeall yo SS 8 (ydlize Cypod) Qyl> clin agi UXy Auobyll lela! os 
Hud 9 BS 9 HSS : blind! Gb b yo SK be ytoine HY [a,b] 


blind! Gb Ad! OVoleoI! plas ALDI ALG Soleo elo i Ao GYolee dotidell Glo 
pg yiiys Ais bb alosl,iI! 5 ball 


64R20 ,65F10 ,65R10 ,45F05 ,47G20 ,45B05 :SWwebyll audlgo aici 


Notations 


R: Set of real numbers. 

C: Set of complex numbers. 

C°{a,b]: The Banach space of continuous functions. 

C'Ja,b]: The Banach space of continuously differentiable functions. 


W??la,b], p € [1,+o0[: The Sobolev space, which vector space of functions that have weak deriva- 
tive. 


H"'{a,b]: The Sobolev space with p = 2. 


LP a,b], p € [1,+00[: The vector space of classes of functions whose exponent power p is integrable in 
the Lebesgue sense. 


T~': The inverse of operator T. 

re(T’): The resolvant set of T. 

R(T, A): The resolvant of operator T. 

sp(T’): The spectrum set of T. 

B(X,Y): The space of linear and bounded operator. 
Ar: The block operator matrix. 

(-,-): Scalar product. 


(-,-): Dual product 
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Introduction 


An integro-differential equation can be defined as an hybrid of an integral equation and an ordinary 
equation or partial differential equation. In the simple sense, it is the relation between the integral and 
the one or more unknown function u derivative’s which appear out side or underside the integral sign. 
Recently, integro-differential equations play an important role, as they are considered the best model 
to express many phenomena in many fields. Its first appearance was a qualitative leap in the field of 
applied mathematics and mathematical modeling. So let’s go back to the early 1900s, when the first 
appearance of this equation was by Volterra, who used it to express the development of individuals over 
a period of time [81]. Volterra may have been the first to use it, but its wide diffusion, development 
and different forms are due to many researchers such as Abel, Lotka, Fredholm, Maltus and Verluslst([?]. 
Later, many mathematicians followed their approach, and many competed either to develop others forms 
of this equation or to devise methods that could help solve it. 


Its first appearance was in biology as mentioned earlier, but this does not prevent us from converting a 
differential equation into an integro-differential equation or from appearing as a solution to a differential 
equation. Its appearance has not been limited to mathematics in all its branches, but has transcended 
many fields. The integro-differential equations are used to model several phenomena [89] such as the 
system of leaking aquifers [29], the earthquake [59]. Also, they are used to define some procedures in 
chemistry or to explain the RL model in electric circuit [55, 54]. The use of these equations are not 
limited to the above domains. It can be used in medicine to model many diseases [28]. The most famous 
of which is malaria [73], the process of coagulation [55], cancer chemotherapy [51], the sugar quantity in 
blood [56] and they are used to model the behaviour of corona virus [35, 85]. 


Throughout history, integro-differential equations have seen different forms, models and types, which 
makes it difficult to solve and to find a single numerical method for different types. So far, many research 
papers have been published with the aim of creating numerical methods for a specific type or improving 
existing methods in order to find numerical solutions closest to the exact solution. They invented and 
built many numerical methods, its main objective is to find the best approximate solution for this type 
of equations, and below we will present some of them: Block boundary value method [90, 11, 17, 86, 50}, 
Jacobi and Gauss-Siedel Method [41], Best approximity point result [14], Operational matrix [37, 64], 
Modified Brensten-Kantorovich operator [12], Hybrid Legendre polynomial [43], Tau numerical method 
[30], Haar wavelet method [42], the variational iterative method [62, 45], Enbedded Pseudo Rung-Kutta 
method [76], ... 


Given the great importance of these equations, many researchers have played an important role in 
studying different models. The most well known integro-differential equation which was studied in dif- 
ferent papers [86, 64, 65, 66] has the following form 

Vu € X, Va € [a,b], Au’ (x) = / K(a,t,u(t)) dt+ f(a,u(x)), u(ao) = uo, (1) 


which is called the Volterra non linear integro-differential equation and X is a Banach space. Another 
equation has a form different to (1) is well studied in [25, 60, 59]. It is given as 


Wwe X Vs ea Fl w= / Kt, @O) 4 FO). (2) 
If x of (4) is fixed in [a,b], we obtain the following non-linear Fredholm integro-differential equation 
b 
Wwe X wes / Kind) al) WF). (3) 
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The linear form of (3) is 


Vu € X, Vax € [a, }], = [ Keo t) dt + f(a). (4) 


We have another equation which related between (2) and (3) is known as Fredholm-Volterra integro- 
differential equation [77] and [73] 


Vu € X, Au(a )= f Kute.t.u(o) af Ko(a,t,u'(t)) dt + f(z). (5) 


The combination of the previous equations allows the birth of a new equation of the following form, to 
which we will give much attention and care in this thesis. 


b 
Vu € X, Au(a =f Ki (a, t)u (ars f Ko(a, t)u'(t) dt + f(x), (6) 
where, in all above equations \ is a complex parameter and f is given function in X. 


Maybe there is no problem in the field that can be modelled in the form of the above equation (6), but 
the mathematical importance of this equation appearance when we use the linearisation method to treat 
(2) and (3) like Newton-Kantorovich method. Therefore, our research can be a reference that facilitates 
the study of this type of equations because through this thesis, we are trying to deal with the different 
forms of this equation. It should be noted that we have searched a little and we find on the other side 
the biggest competitor of these equations, which are the differential equations. We find that any problem 


of the following form 
Au!" (a) + a(a)u'(x) + b(ax)u(a) = g(x), 
u(a) = a4, u(b) = fr (7) 
u'(a) = ag, u '(b) = Ba, 


can be written in the previous form (6). 


The use of differential equations is clearly shown in physical phenomena of solids in free vibration, in 
electricity through a circuit in mechanics through Newton’s equation and many examples shown below 


Simple Mass-Spring Systems a 


Spring: 
Elastic beam 


Masses: 
Masses of the 
bridge structure 


Figure 2: Physical demonstration of a simple mass-spring system: A mass attached to an elastic cable 
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Motion 


Disturbance causes 
vibration of the mass 


Figure 3: Physical demonstration of a simple mass-spring system in free vibration 


The last examples have the next form: 
u’(t) + Lave = 
m 
where, u(t) signified the instantaneous position of the mass, m is a mass and & is a spring constant. 
Damper- Spring: 
a dashpot: k —_ 


Cc 


Dashpot for 


Figure 4: Physical model of damped vibration 


This example has this form 
c 


k 
A / 
u(t) + me (t) + mult) = 0. 
Therefore, it can be written in the form of the above equation (6), and by using the clear methods of this 
equation, you can get better numerical solutions than the existing numerical solutions. Here, we have 
clarified the importance of the equation (6) and we are going to study numerically and analytically. Now, 
we would like to review and discuss the focus of our thesis. 


Our thesis contains four essential chapters. In the first chapter, we give some concepts and definitions 
that we are going to use in the thesis. In the second chapter, we study and to analyse the existence 
and uniqueness of the solution in the Banach space X = C1[a, b]. We construct three different methods 
based on the Nystrém, collocation and Kantorovich methods to approximate the solution. We introduce 
theorems that show the convergence of numerical solutions. At the end of the chapter, we compare the 
numerical results of the three approximated solutions and develop a new iterative system since the system 
obtained by the three methods is of very large size. 


In the third chapter, we study analytically and numerically the equation in the Banach space X = 
Ca, b| but this time the two kernels of our equation are weakly singular. So, we find sufficient conditions 
to prove the existence and uniqueness of the solution. For the numerical approach, we construct two 
different methods. The first one is based on the product integration method. The second is based on 
b-spline collocation methods. The difference between the two is in the size of the system and the speed 
of convergence. While the first method gives us a system of four blocks, the second method gives us a 
system of only one block. The consistency of the numerical methods is well explained and the comparison 
between the two methods is well illustrated in the numerical example. 


Because, we are interested in the smallest details of our equation, here we study the uniqueness and 
presence of the solution in the most complex Banach spaces: Sobolev spaces. Especially since they are 
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considered to be the best spaces that well express many real-life phenomena. We study the uniqueness 
of the solution and its existence in the space W!’[a, b], but the numerical methods developed one only 
W!la,b] and H+[a,b]. Our numerical solution will be obtained by constructing two methods based on 
each of Kantrorovich and Galerkin. After having obtained our apical solution, we study its convergence 
through several theorems which explain and prove it in the sense of the norm of two spaces W1![a, b] 
and H'{a, 6]. All this will be expressed in the last chapter . 
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1.4.2 Approximation Methods For Fredholm Integral Equation with Weakly-Singular 


Kernel... . 


In this chapter, we shall begin by introducing the essential basic notions that we use during this 
thesis. We start by recalling some notions about the functional spaces. Then, we present some notions 
and theorems on the linear and bounded operators which helps us to introduce the classical process of 
integral equations defined on a bounded interval. We focus on the linear Fredholm integral equations 
and we propose some approximate methods applied to this type of equations. The presentation of the 
Fredholm integral equations allow us to well introduce the integro-differential equations concept. Note 
that, all theorems and properties will be presented without proofs. 


1.1 Operator Concept 


Let X and Y be Banach spaces with the norms ||.||x and |].||y respectively. 


1.1.1 Linear and Bounded Operator 


A function T which maps X into Y is called a linear operator [23], if 


Vu,we Xanda,y€C(orR), Tl(av+ yw) =aT(v) + 7T(w), 


and bounded if there exists a positive constant (, such that 


Wo € X, ||Tolly < Allellx- 


The norm ||.|| of the operator T is defined by 


[PI] = 


[Polly _ 


sup ||Tully = sup ||Tolly. 


Ilelix%o [ellx jolie <a lolx =1 
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1.1. Operator Concept 


Theorem 1.1.1. /2/ Let T : (X,]||.||x) > (¥Y|.||y) be @ linear operator. Then, the following assertions 
are equivalent: 


1. T is continuous on all X, 
2. T is continuous in 0, 


8. T is bounded. 
Theorem 1.1.2. The linear space B(X,Y) of bounded linear operators from X into Y is a Banach space. 


Notation 1.1.1. 
1. We denoted by B(X) the Banach space of linear and bounded operators from X into itself. 


2. The operator that sends each v € X to 0 will be denoted by O and the identity operator will be 
denoted by I. 


Definition 1.1.1. /2] Let (Tn)nen € B(X,Y) an approximations sequence of T. (Tn)nen converges 
pointwise (or simply) to T, if only if 
lim ||(Tn—T)vlly =0, Vue xX. 


n—-+oo 


Definition 1.1.2. /2] (Tn) nen ts called a uniform approximations sequence of T € B(X,Y), if and only 


if 
lim ||Z, —T|| =0. 


n—->+oco 


Theorem 1.1.3. (Banach Steinhaus)/10/] 
Let (Tn)nen be a family of linear and continuous operators in X from Y. Assume that 


sup ||T,vul|ly < oo, Vu € X. 
neN 
Then, 
sup ||Tin|| < 00. 
nen 
In other words, there is a positive constant 8 such that 


Vu € X, ||Thully < Bllvllx, Vn EN. 


Corollary 1.1.1. //0] Let (Tn)nen be a sequence of linear and bounded operators from X to Y, such 
that T;, converges pointwise to T. Then, 


1. sup ||Tn|| < co, 
nEeN 


2.T € B(X,Y), 


3. ||T|| < lim inf ||T;|. 
n—+o0o 


Remarque 1.1.1. 


e The second result of the last corollary shows that (Tn)nen converges to T uniformly on compact 
sets. 


e Inthe American literature the Banach Steinhaus theorem’s is referred to as the Principle Uni- 
form Boundess, which expresses the result: We deduce a uniform estimate from point estimate. 
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Operator Concept 


1.1.2 Compact Operator 


In this section, we present the notion of a compact operator with some theorems concerning this type 
of operators. 


Definition 1.1.3. /3// Let T © B(X,Y). It is called compact if it maps each bounded set in X into a 
relatively compact set in Y. 


Lemma 1.1.1. [/5] Let (fn)nen be an equicontinuous sequence in C°la,b] and f € C°[a, b] such that for 
each x € [a,b], |fn(a) — f(x)| > 0 as n + +00. Then, ||fn— fll 30 as n— +00. 


Theorem 1.1.4. /84] 


1. An operator T € B(X,Y) is a compact operator if and only if for every bounded sequence (Un) nen € 
X, the sequence (TUn)nen has a convergent subsequence in Y. 


2. (Arezla-Ascoli theorem) 
Let D be a compact set and M a subset of C°(D) with ||.||.. M is pre-compact if it is 


(a) Uniformly bounded, i.e sup || flo <0. 
feM 


(b) Equicontinuous, i.e; Ve > 0,2 € D,A6 with, 
lf(z)-—fly)| <e, VfeEM, andy € D, with |x — y| <0. 
3. Compact linear operators are bounded. 


4. Linear combinations of compact linear compact operators are compact. 
Definition 1.1.4. //0] An operator T € B(X,Y) is called the finite rank if dim (Ran(T)) < +00. 


Theorem 1.1.5. //0] 
1. Let T: X —Y is a bounded linear operator of finite rank. Then, T ts compact. 


2. If (In)nen € B(X) is a sequence of compact operators and \|T, —T\|| > 0 as n > oo. Then, T is 
compact. 


1.1.3. Spectral Theory 


The spectral theory plays an important role in our study. So, in this part we give some definitions and 
theorems concerning this theory. 


Definition 1.1.5. /2/ 


e The kernel ( or null space ) of a linear operator T is 
Ker(T) = {v EX:Tv= of. 
e The range of T is the image of X under T: 


Ran(T) =T(X) = {we ¥, w= Tov, we xh, 


Proposition 1.1.1. /2/ 
e The operator T € B(X,Y) is injective if Ker(T) = {0}. 
e The operator T € B(X,Y) is surjective if Ran(T) =Y. 


e The operator is bijective if it is injective and surjective. 
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Definition 1.1.6. /2/ The operator T € B(X,Y) is inversible, if T is bijective and T~' is bounded. 


Definition 1.1.7. [2] Let T © B(X) The resolvent set of T is defined by 
re(T) ={AEC: AI —T is bijective }, 


and for \ € re(T), R(T, ) = (AI —T)~! is called the resolvent operator of T at X. 


Definition 1.1.8. /2] Let T € B(X). The spectrum of T is the set 
sp(T) = C\ re(T),= {r €C: AI—T has no inverse \ 


and an element of sp(T) will be called a spectral value of T. But, A € sp(T) is called an eigenvalue, if 
the equation Av = Tv has a non null solution (eigenfunction ) v € X. 


Theorem 1.1.6. (Neumann Expansion)[13] If T is a bounded linear operator from X into X and 
A €C such that ||T|| < |A|. Then, AL —T has an inverse in B(X,Y) given by the following uniformly 
convergent series 


+00 T 
RUT.) = (AI- T+ =)" ee 
k=0 


and 1 
Ota) ||S5— aa 
|A| — ||Z| 


Let T € B(X). If \ is a given complex paramater. Then, for every w € X, a solution v € X of the 
linear equation 
AI-T)v=u, 


will be uniquely determinate by w if and only if \ € re(T) and v is given by 


v= R(T, A)w. 


1.1.4 Collectively Compact Operator Approximation 


The notion of collectively compact operators is important to show the convergence of numerical solu- 
tions. For that, we introduce it in this section. 


Definition 1.1.9. /75] The sequence (Tn)nen € B(X) ts called a collectively compact approximation 
sequence of T € B(X) if and only if (Tn)nen converges pointwise to T and there is exists no > 0, such 
that 


U {(a=Tye ve X, |lv|lx < i} 


n>no 


is a relatively compact subset. 


Theorem 1.1.7. [/8] Suppose that (Tn)nen is a collectively compact sequence in B(X). Then, the 
following hold: 


1. Each T,, is compact, 


2. If T, > T pointwise. Then, T is compact. 


Definition 1.1.10. /2) The sequence (Tn)nen € B(X) is v-convergences to T € B(X) if ||Tn|| < +00 
and 


(Tn — T)T|| =0, (In — T)Ty|| = 0. 


lim || lim || 
n—->+oo n—->+oo 
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Operator Concept 
Theorem 1.1.8. [75] Let (Tn)nen € B(X,Y) be a sequence of operators pointwise convergent to T € 
B(X,Y) and A is a compact operator. Then, 

lim ||(Z; — T)Al| = 0. 


n—>+oo 


Proof. Since A is a compact operator, the following set 


ee {Ae llollx < if, 


is relatively compact in the Banach space Y. By the Banach-Steinhaus theorem 1.1.3 7;, converges 
uniformly in Y and 


lim ||(T, -—T)A||= lim sup ||(Tn —T)Av|ly = jim sup ||(T, — T)w||y = 0. 
a K 


n—+00 n—- +00 llvl|x<1 SO wE 


1.1.5 Block Operator matrix 


A block operator matrix is a matrix, where its elements are linear and bounded operators. Every 
linear and bounded operator can be written as a block operator matrix if the space in which the operator 
defined is decomposed in two or more components. 


. N 
Definition 1.1.11. Xj = |] X; is the Banach product space, equipped with the following norm 
i=0 


N 
VV = (v1, 02,---0n) € Xn, [IV1lx = dolleillx:- 


i=0 


Definition 1.1.12. Ar is a block operator matrix with the following representation 
Ar : xe '—> Xn 
N N N 
V r> ArV = (Sos), DoFai05,--- Tse Jal <jcn, 
j=0 j=0 j=0 


where, {Tijhi<ij<n are linear and bounded operator. 


Definition 1.1.13. B(Xn) is the Banach space of linear and bounded operators defined in the product 
space X into X, with the norm 


V Ar € B(X),|I|Ar|l| = sup ||ArVIlx,- 


Vile yaa 


Definition 1.1.14. Let I be the identity operator and O the null operator of the Banach space X. We 
define the block identity operator as In : Xn —> Xn 


Ir oO oO O 

Oo. ft £ O 
Iy = 

O 0 . O F 


Proposition 1.1.1. The norm of the block operator matrix Ar is given by 


2 
Aci = may) I Zol 
q= 
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1<p<N 


Proof. We have |||Ar|||= sup ||ArU||x,, such that, 
U yal 
N ON 
ATU ley = DOM OTeellx, 
q=1 p=l 
N N 
= SS llZevtiell xs 
q=lp=1 
N N 
& Sy ell Ilupllx; 
q=1p=1 
N N 
S DillZerlldollupllx, 
q= p=1 
N 
< max S7l|Tupll lx: 
q=1 


This gives 


N 
es S- 
I|Ar||| < Rae ‘Feel 
q= 


N N 
Without loss of generality, assume that max > ||Tqp|| = >> ||Tqill- 
ls<psNq=1 q=1 


Let Up = (1,0,...0). We have 


N N 
\|ArUoll = So ||Taall = Be, 2 lITenl 
= 


q=1 


Then, 


N 
HA] = nar) Zl 
q= 


1.2 Cubic Spline Interpolation 


Cubic spline interpolation is considered as a problem of drawing a smooth curve through a n+1 points. 
So, we can use a polynomial of degree n > 1 passing through all these points, but this is not available if 
n is large enough. For this reason, we use different polynomials in different regions. Then, the obtained 
interpolation function is piecewise polynomial. In addition, these approximations may be continuous but 
in general the derivatives can be non-continuous. 


1.2.1 Smooth Cubic Spline 


The essential purpose of constructing spline functions is to find an approximate derivable functions. 
So, let [a,b] be an interval of R, let the given function f(x,;) for i = 0,1,...n, such that {a;}"_) are 
defined as a = a < 4 <-+-+: < 2, and let g be an interpolation function of f on the interval [a,b]. We 
divide the interval [a, b] into n subintervals I; = [x;-1, v;] for 0 <i <n. The first simple approximation g 
can be fined by employing the linear interpolating in each subintervals J;. In this case, f is approximate 
by a straight line. This approximation g is refried to linear spline which is not be smooth. To fined a 
smooth approximation, we can use a polynomial of higher degree over each of the subintervals [;, such 
as the piecewise cubic interpolation. 


In cubic spline interpolation, we try to find an approximation S for any given function f, verifying the 
following conditions 


1. S(x) € C?[a, bf, 
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2. S(a;) = f(a), fori =0,1...n 


3. On each interval [;, S(x) is a polynomial of degree 3: 


Si(z), oS U< 4K, 

Siz), trex x, 
S(x) = . 

Sn(x), Ci ES ey 


where each S; has the following form 
Si(x) = a4 + ja + yx? + dja?, i=0,1...n. 
To determine the cubic spline S(x), we need to calculate the coefficients a;, b;,c;,d;, by assuming that 


i-1)9 = f(yw-1i),i=1,...n 
Si(z) = f(a), i=1,...n 
) 

) 


II 
RR 
+ 
iG 
ny 

8 
— 

© 

| 
= 

3 

| 
i 


We have 4n — 2 conditions. But we must calculate 4n coefficients of the spline polynomial. Then, we 
need to add two conditions and there are three types of these conditions like 


1. Si (x0) = f (xo) and Sylkn) = f' (Zn), 
2. Si/(ao) = f"’ (xo) and S’ (an) = f" (an), 
3. Si(xo) = Sn(xo), $1 (a0) = S’n(an) and $”1(20) = Sn(an). 


There are many methods to calculate the coefficients of the polynomial spline. But, we choose the fol- 
lowing explicit form of cubic spline based on the moments M; = S$” (a;) 


Ee 3 ,. -p2 a 
Si(a) = Mi-1 6h Gh | (se) wat \(& *) 


si(or) = Bas + Ban + fed fle), 
Si(az) = aM, * Mey Move F(ei) 


The continuity of S/(a) at 2; yields 


Vika ANE xs = 6( Hee Ae Ae) fe), i=1,2,...n-1. 


h2 h2 


But, for many applications, it is more convenient to work with the slopes m; = S’(#;) rather than the 
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moments M;. Here, we present another interpolation in each segment [x;_1, x;] as 


Si(z) = mi-a (x: = 2) {a = i=) Mi ta mo) (xi — 2) 

& seo ace tina) +A) pcp, (@ = Bi-1) Be = 2) +h) 
Sie) = mia (ai — = i = 3a) ae (a — sees vi-1 — 32) 

_ fi) P(@i-1) 6, acs 


T h3 


The limit values of second derivative at x; are 


= 4 Xi) — bi 
Si(z7) = oe es ofl ) a ) 
—4 2 f(vi41 — f (xi) 
+) _ = 
Sia) = lL eae + nT 
For finding the results, we require solving this system 
ati i=1,94.40=— 1 


h 3 


1.2.2 B-spline functions 


In many applications it is necessary to express the approximations functions as a linear combination 


of approximate basis functions.This is case for differential or integral equations. 


A simple presentation of spline in terms of independent basis functions can be given as: 


m—n—1 


Yn >1, Va € [a, 6], > a; By n(2). 
i=0 


where the m — 1 b-splines of degree n are fined by the following recurrence formula: 


L- 2; 
Vane a, bl, Bi n(2) = ———*_ B; n-1(2) + 
Titn — Vv Titnt1 — Viti 


Xvi — 
ee Begg ie) 


with 
’ if UCL XY ’ 
Bi o(x) = { 0 else. 7 


L— Li-1 A 
; ifaj1 <2 < 4%, 
Uj U-1 
: = Liste , 
Bya(2) kad sel ifa; <2 < 241, 
Ti41 — 
0, else. 


It is clear that the latter (1.1) is a linear piecewice function. 
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For n = 2, we have 


(a = rae 


U1 U< 2; 
(ti41 = Li-1) (Xi = Li-1) ‘ . 


d 


(ote 4) (eg — ar] (tit2 — 2i)(@ — i) me < eee 
iS i+1 
(Bigg = Bi) (G44 a te) (Xi42 _ Hoe Cea = a) : ial 
Bia(x) = ; (12) 
(zi42 — 2) 
: < < 
(win = 24) (Xi42 = isi)’ TMi SUS LH42, 


0, else. 
which is called the quadrature b-spline basis. 
For n = 3, we get 


(a = a1)? 


Tp SE << G 
(Xi42 = Li-1)(Li41 = Li-1) (Xj = Xi-1) : 7 


b 


(x — %-1)"(wi41 — 2) _ (© — ~1)(®i+2 — 2i)(@ — i) 
(Xi42 = Li-1)(Vi41 = Li-1) (Lig a xi) : (Ui42 — Li) (i42 _ Xi) (i441 — xi) 


(x — %i~-1)(ti42 — 2)” 


Tee <i 
| (Li42 — Li-1) (Via _ Li) (Li42 = aaa) ‘ “ 


(vi43 — x)(x — xj)? (te. )\(tee— 97)" 
(ri+3 - ti) (Bare = Ge aay _ %;) (ri+2 a ti1) (Zi = Cty = Tigi) 


(vi43 — ©)(@ — 2i)(Zi42 — 2) 
(xi43 = He) (Bete = Ba) Bis = 2441) , 


+ DiS ES Pe; 


(ti3— z)? 
(@ir3 — 2) (rs — By (tg — aye)’ 


0, else. 
(1.3) is a cubic b-spline basis. 


1.3 Continuity Modules 


In the thesis, we need to use different formulas of continuity modulus. In this section, we present a list 
of these formulas. 
We start by defining the continuity module in X = C°{a, b] 


VR >0,VVEX, Ko 0(9,h) = sup |v(x) — v(y)|. 
|z—y|<h 


For all functions define on the square [a, b]?: 


VreE [a,b], Vh> 0,Vg € C° ([a, 8]?, R), Ko0,0(g, 2) (2) = sup lo(x, 41) _ g(x, y2)|, 


|y1—y2|<o 
and the continuity module k.o,1 as 
0 
V@eE [a, b), Vh > 0, oe € C°([a, b]?,R), Roo, 1 (9; h)(x) = Ro0,0(9; h)(x) = K0,0(Org, h)(x), 


where Ozg is the partial derivative of g respect to x. 


Study of linear integral differential equations of the second kind: Analytical and numerical approach 17 


1.4. Fredholm Integral Equations 


We present the continuity module in the Banach space Xz = X x X as: 
Vh>0,VV= (v1, v2) € Xo, Koo 1(V, h) = Koo,0(U1, h) + Koo,0(U2, h). 


We introduce the continuity modulus wo(h,.) in L*[a, b] by 


Vg € L*[a,b], K1,0(h, v) = ow | |g(a + y) — g(y)| dy, 
x€E[0,h 


and the continuity modulus «1,1(h,.) in W*[a, b] by 
Vg S W""la, b], K1,1(h, g) = wo(h, v) =r wo(h, 9g’). 


We give the generalization continuity modulus for any continuous function G(z,t) € L'({a, b]?, R) on the 
square [a, b]?: for any t € [a, | 


F10(h,@)(é) = sup [ IG + y,t) — G(y,8)| dy. 
x€([0,h] 


G 1 2 
Fe (tet) € L'([a,0?,R) 


We introduce the generalization continuity modulus for 


Vt € [a,b] K11(h, G)(t) = K1,0(h, G)(t) + K1,0(h, OG) (t), 


where 0,G = (6 .), and we define 
xv 


[1,102 G)lloo = mane |fia,1(h, G)()|- 


1.4 Fredholm Integral Equations 


Many problems and equations in physics can be transformed into Fredholm integral equations like the 
initial and boundary value problems. Every Fredholm integral equation contains a function obtained by 
the unknown function u by integration and of the form i K(a,t)u(t) dt, where K is called the kernel 
and is assumed known. In general, the Fredholm integral equation contains the known function called 
the free term. Usually, they are complex-valued functions of real variables. 


This type of integral equations has three kinds and each kind refers to the localisation of the unknown 
function u. First kind has the unknown function present under integral sign only 


b 
x) = K(a,t)u(t) dt 


The second and third kind have the unknown function outside the integral sign. Thus, the following 
equation is of the second kind 


n= [ Ke.ou0 ) dt + f(a), (1.4) 


where A is a numerically parameter generally complex. In particular applications \ is composed of physical 
quantities. 
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In this part, we set X = C° a,b] with the norm 


Vu € X,||v||x = max |v(2)|. 
a<a<b 


We define the linear integral operator T by 


T: X —>X 


b 
u > Tu(x) = / K(a,t)u(t) dt, a<a<b. 


Then, the integral equation of the second kind (1.4) is rewritten as 


(Al —T)u = f. 
Recalling that 
[Z|] = sup ||Tullx, 
|u|] <1 


the norm of operator T is given as 
b 
TI|= max, [1K (2,1) at, 
acxcb Ja 
because for each u € X with |]u||x <1 
b 
|Tu(x)| < max / |K(a,t)| dt, Va € [a, bd], 
acacdb Ja 
and so 
b 
\|Z || < max i |K(a,t)| dt. 
acrcb Ja 
Since K is continuous, there exists x9 € [a,b] such that 
b b 
[1K o.nlae= max, | |K(@.0) at 


We set 


sup _||Tul|x = ||Tuollx- 
[lull <1 


Then, for all ¢ > 0, we choose up € X by 


= K (ao, t) 
~ |K(xo,t)| +e 


uo(t) t € [a, 6). 


It is clear that ||uo||x <1. So that 


b 2 b 2 2. 
| F< (xo, t)| / | (xo, t)| —e 

T — + dt > —_—__—__ dt. 
|Puo(eo)| a eecies: = J, [Kot] +e 


Hence, 


b 
ini f |K (o,t)| at — e(b — a). 
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Since this holds for all ¢ > 0, we have 
b b 
> = : 
I> f ie o.0|de= max, [x 2,0) a 
We say that if ||T'|| < |A] the Fredholm integral equation (1.5) has a unique solution in X. 


1.4.1 Approximation Methods For Fredholm Integral Equation with Contin- 
uous Kernel 


In all most common situation where an exact solution to a problem cannot be found directly but 
its existence and uniqueness are assured as in the previous case, the approximation techniques become 
important. In what follows, we discus some methods which can deal with successfully. 


Projection Approximation Methods 


We are searching for an approximation solution of the integral Fredholm equation (1.5) based on the 
projection method principle. We need to define the sequence of finite rank (P,)n>1 € B(X). We de- 
note by X,, the sequence of finite dimensional subspaces of the Banach space X which have the basis 
{e;}7_5, by X* the dual subspaces of X,, with the basis {e; }'_) and by (.,.) the duality bracket, such that 


fl if i=j; 
(ene) = { 0 if iAj. 
The sequence of finite rank (P,)n>1 € B(X) is given as 
Pr: Xn 7X 


n 


u + Pru => (u, ef es. (1.6) 
i=0 


This approximation operator is used to construct several projection methods like Galerkin and Kan- 
torovich methods which will take a large part in our study. 


Lemma 1.4.1. Let T € B(X) and (Pr)n>1 € B(X) be a sequence of finite rank presented by (1.6). If 
(Pr)n>1 converges pointwise to the identity operator and T is a compact, then 


jim || — Pr) Tl] = 0, 


Proof. See theorem 1.1.8 proof. 


Theorem 1.4.1. Let T € B(X). Assume that (AI — T)~! exists, is bounded and 
jim || — Pr)T|| = 0. 


Then, for alln >1,(AI — P,T)~! exists from X to X. Moreover, it is uniformly bounded: 


sup ||(AI — T)~*|| < 00. 


n>no 


Proof. For n > 1, we have 


M-P,T = (Al-T)+(7-B,1), 
(AI — PaT)[T + (I — P,T)71(f — PpT)]. 
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Now fix no € N such that 


1 
—— By ee! ies Sr = eee 
a ee Oren 


Then, 
WAL — TT) (L — PaT)|| < |JAL-—T)~*|| [If — PaT)|| <1, 


which gives us by Neumann’s theorem 1.1.6, the existence of these inverse 


[r+ Or-T) (rT -—P,7r)-}, 


which is uniformly bounded 


=1 -1 1 
I+ OF-T'0 - PANW S GAT 


All of this proved that the inverse (AI — P,,A)~+ exists and can be written as: 


Ql-—P,7T) + =[F+Ql-T) 17-27) Or-T), 


with 


AL=7)7"| 


M=P,TY “|= = 
i NS Ty ell F— 2) 


Ch. 


Galerkin Method 


The word Galerkin is used in different ways by researchers in books and research papers to define 
Fredholm’s integral equation approximation solution. For example, in Atkinson [5] and Nair [48] books, 
Galerkin method is a name of a numerical treatment used to find an approximation sequence (Un) n> 1 
solution of 


(AI — PpT)un = Ppf. (1.7) 


However, in Ahues in [2], the principle of Galerkin method is employed to determine the sequence (Un)n>1 
solution of 


(AI — PaxTPy)tin = f. (1.8) 


We focus on the definition (1.7), for this reason we give the following theorem 


Theorem 1.4.2. Let (Pn)n>1 be a sequence of finite rank projection operator defined by (1.6). The 
sequence (Un)n>1 given by 


n 
Un = s aei, 
i=0 


is a solution of (1.7) if and only if we get the following system 


day = Yoai(Tei,e}) + (fe9), 9 = 0,1... oe 
i=0 


Study of linear integral differential equations of the second kind: Analytical and numerical approach 21 


1.4. Fredholm Integral Equations 


Proof. Assume that 


Un = S aie: 
i=0 
From the definition of (Pp )n>1 in (1.6), the equation (1.5) is equivalent to the following 


Min = SI Tun, €; )e; + f,e;) 
1=0 
Multiplying by ej, we get 
A(Un; ej) = (Tun, ej) + (f, 5). 


From the equation (1.10), we get 


n 


A(Sravese}) = “(Pawn ‘je + he ej Jee, 
i=0 = 


i= 


which is equivalent to 


‘Yea (€;,€ 9) = Yoo (Te;, e% €;) + (f, e;). 


Finally, 


Aaj = S “ai (Tei, ef) 4 fi) 


Theorem 1.4.3. Let un and u be solutions of (1.7) and (1.5) respectively. Then, we obtain the two-side 
error estimate 


al 
[AT — PTI 


||u — Prullx < |lu—tallx < JA] AL — PrT)~*|| [lu - Paull, 


and if (Pn)n>1 converges pointwise to the identity operator, then 


lim ||u — up||x = 0. 
n—->co 


Proof. Applying the projection operator P,, on the equation (1.5), we obtain 


AP,u — PyTu+ Au — Au= Pyrf, 


which equivalent to 


(AI — P,T)u = Prf + A(u— Pru). 


Substracting (AI — P,T)u = P,f to get 


(AI — P,T)(u — un) = A(u — Pru). (1.10) 


Study of linear integral differential equations of the second kind: Analytical and numerical approach 22 


1.4. Fredholm Integral Equations 


Then, 
Iu — unllx < [Al [AL — PrT)~*I| |lu — Prullx- (1.11) 
From (1.10), we have 
Al [lu — Prullx = |[(AL— PpT)]| [lu — unl lx- (1.12) 
Moreover, 
WAL—PaT)|| < |AL- T+ - Pa) TIL, 
< ||AL-T||+ én. (1.13) 
Substituting (1.13) in (1.12), we obtain 
I le Paull < (IA — TI] +60) lle — unl (1.14) 


Hence, 


|A| 
AL— Tl) +én 


Iu — Prul|x < |]u— unllx- 


From (1.11) and (1.14), we get the result. Moreover, we have 


lim |Ju — P,ul|x = 0. 
n+ co 


This implies that u, converges to u in the sense of the norm of Banach space X. 


We shall consider two special case of Galerkin method, depending on the framework space. The first one 
is collocation method if X = C°[a, b] and the second one is the orthogonal Galerkin method if applying 
this method in Hilbert space X = L?{a, b]. 

(i) Collocation Method 


We take X = O° [a,b], and we define A,,, for all n > 1 the uniform discretization of interval [a, b] as 


An = fn 1a = ty <a << ty <n, h = 2541-25, o<jnk (1.15) 


Now, we return to the essential condition of the basis’s choice (1.6). For this reason, {e;}"_) are 
the following haat functions 
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1+ aes iad yt |; 
ei(x) = 
0, 
,a22 x E [x0, x1], 
_ h 
eo(x) = 
0, 
i LE [Zn=1;2nl, 
_ h 
n(x) = 
0, 


with the duality bracket is defined as 
(u,e;)=u(a;), ue xX,i=0,1,...n. 


Then, the sequence {P,,}n>1 verifies the following interpolation condition: 


P,u(2;) = u(a;), +=0,1,...n. 


It is clear that {Pn}n>1 is an interpolation projection based on nodes 29, 21,...,2n.- 
The system (1.7), is equivalent to 


AX =MX+F, 


where X = (a9,Q1,---Qn)' € R"*!, F = (f(a), f(a1),---, f(a@n))* € R™t! and M is a matrix of 
size (n+ 1) x (n+1), such that its elements are given by 


1 xy Li4. 
(Tei,e5) = ; | K (xj, t)(t — 4-1) a+ | K (x;,t)(wi+1 — t) dt, O0<j<nl<i<n-l, 
Ti-1 Xi 
1 f* ; 
(Teo, €;) _ xf K(2;,t)(a1 — t) dt, O<j<n, 
xo 
(Ten, e}) = i K(a,;,t)(t — tn-1) dt, O0<j<n. 
Ln—1 


(ii) Orthogonal Galerkin Method 


Assuming that X = L?[a,b] with an inner product (.,.). We recall that the duality bracket is the 
inner product in this Hilbert space. 
For each n > 1, let X,, be a finite dimensional subspace of X. Let {e;}%_, be a basis of X,, checks 
the condition (1.6), that means it is an orthonormal basis with respect to the inner product (.,.). 
Then, the sequence (P,,),>1 is an orthogonal projection. Now, it has the above presentation 


Vue X, P,au= S"(u, ei)ei, 
i=0 
and the system (1.10) has this form 
day = Dros (Tes,e7) + (Frey), J=0,1,...m, 10) 
i=0 
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System (1.16) can be written as 
\X =MX+F 


where X = (a9, 0,..-Q»)* is a vector of R"t', M is a matrix with size (n+ 1) x (n+1) such that 
each element mj; are given by 


b pb 
m= ff K(a,tes(teda) dt dr 0<ij <n, 
and each element F;,i = 0,1...n of the vector F is given by 
b 
F; =| f(ax)e;(x) da. 


Kantorovich Method 


In Galerkin method the convergence of the approximate solution to the exact solution is guaranteed by 
the convergence P,,f — f. To avoid this situation which is considered as an inconvenience, we propose the 
Kantrorovich method where this approximation procedure was first studied by Schock. For this purpose, 
we assume that f € Ran(T). 


Recalling the Fredholm integral equation 


Yue xX, Au=Tut+f, 


and applying the linear integral operator T to both sides of equation, we obtain 


Av=Tv+f, 
with v = Tu. Now, we call u* a Kantorovich approximation of u and it is defined by 


1 
uf = 5 (08 +s), 


G is a Galekin approximation defined in the last section and verifying this approximation equation 


where, v,, 


hw? = P,, ToS + PT f. 


To illustrate the relation between two methods, we observe 


Auk = wo +f, 
1 
= 5 (rare? +P,T ') noe 


1 
= 5 PoTon +f) +f 
Recall that v¢ = Tu. Then, the Kantorovich approximation equation has the following form: 


Mu® = PTuk tf, (1.17) 
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Theorem 1.4.4. Let uX be an approximation solution and u be an exact solution. Then, we have 


ME — PrP) Z 
IAT—TI+0-PaT 


Iu — wt Ix SAT — PrT)"*|I ||P — PaT)|I- 


Moreover, assuming that P,, converges pointwise to the identity operator. We get 


F _, K = 
Jim |lu— Un lx = 0. 


Proof. For n large, we have 


OF-P,T)\(u-u*) = (T-P,T)u, (1.18) 


(QI-T)+(P,T —T)|(u-u*) = (T-P,T)u. (1.19) 
From (1.18), we have 

leur lle < IAP- PP) IE = PPL; 
and from (1.19), it follows 


MT — Pa) - 
IAr=T+1@-P.7) = 


Iu — un Ix. 


Then, we obtain the result. 


Nystr6m Method 


In this part, we describe another approximation method called the Nystrém method which is highly 
applicable for the numerical solution of the Fredholm integral equation of the second kind (1.4) in the 
Banach space X = C°/a, b. 


We recall A, for all n > 1 the uniform discrtization of interval [a,b] (1.15). This method is based on 
the next numerical integration scheme 


b n 
7 g(x) dx = S \wig(wi), (1.20) 
e i=0 


where {w;}%_p9 are called by the weights under the condition 


su W;| < +00. 
sp il 


The Nystré6m approximation u, is the solution of the following numerical equation 


Va € [a, }], Min (& = Yok (@, Uj )Un (vs) + f(x). (1.21) 


We select the collocation points z;, 7 = 0,1,...,n to obtain the solution at this knots. We get the 
following linear algebraic system 


Atin(X;) = Yok (xj, 04 )Un(ai) + f(z;), j=0,1,...n 
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We return to the equation (1.21) which has this simple formula 
Va € [a,b], Mn(£) = Thtin(x) + f(x), 


with T;, is linear and bounded numerical operator of finite rank. It is defined by 


Despite the wide application of Nystrém’s method on integral equations due to its simplicity, but only 
a number of researchers who developed the error analysis of this numerical treatment during the 1940s- 
1970s year. At the beginning, the aim of the researchers was to demonstrate the stability and convergence 
of the method. Later, the goal changed to the creation of an error analysis in a more abstract framework 
that contains all possible situations and what we mean by them the error domain, the convergence rate 
and the functional space. The concept that generalises all possible cases is the concept of collectively 
compact approximation where the Nystrém approximation operator can be verified by. 
Now, we present the most important theorem that helped for the creation collectively compact operator 
concept approximation. 


The next theorem shows the pointwise convergence of (T;,)n>1 (1.22) to T (1.5). 


Theorem 1.4.5. /5, 48] Let (In)n>1 be the Nystrom approximation of the integral operator T. Then, 


1 U {Ta we X, |lullx < i} is equicontinuous, 


n=1 


2. (Tn)n>1 is pointwise approximation of T. 


Proof. see [5, 48). 


Corollary 1.4.1. Let T be a compact operator and T,, defined by (1.5) and (1.22) respectively. Then, 
lim ||((Z—T;,,)T|| =0 and lim |\(T — T,,)T|| = 0. 
noo noo 


Proof. By theorem.1.4.5 T), is a pointwise convergence of the compact operator T and the set 


M= {Ta ue X, ||] | x Si, nen}, 


is equicontinuous. Hence, by theorem 1.1.8, we obtain the result. 


1.4.2 Approximation Methods For Fredholm Integral Equation with Weakly- 
Singular Kernel 


So far in this thesis, we have been interested in the numerical methods to solve Fredholm integral 
equations under the assumption that the kernel is a smooth function. Among cases the domain kernel’s 
definition is infinite or the kernel has a singularity with its domain of definition. The integral equation is 
said to be singular. 

Let X = C°[a, b], the Fredholm integral equation with weakly singular kernel has the following form 


Vu € X, Au(x) = [ole —t)h(a, thu(t) dt+ f(r), a<au<bd, (1.22) 


a 
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where u is an unknown function. We assume that f € X, H(a,t) € C°({a,b]?,R) and the function p 
verifies the following hypothesis 


p€ L1(0,b—a), 
(99) 
p(0) = 0. 


In general, p has the following form 


1. Logarithmic form 
p(\ — t]) = log |a — ¢]. 


2. Algebric form 
p(ija—tl) =|a—-t*, O<a<l. 
Product Trapezoidal Rule 


Let n > 1, we recall A,, the disceritization of the interval [a, b] defined by (1.15). 
Vue X and x;_1 <t < x;, we define 


[H(a,t)u(t)|n = - (a, — t)H (x, v1) u(a;_1) + (t — a1) (a, 2;)u(a;)}, V +=0,1,...0. 


Let T;, be a numerical approximation of the integral 


Vu € X, Anu(x) = Sow; (x) H(a, xi)u(ai), Vx € [a,b], 
1=0 


with the weights {w;}'_9 defined by 


wo(a) = 7 f plle—a)ler—0) at 


(0) 


1 f* 
5 [eke = tlt = stn) a 


ri 


1 Citi 1 Li 
wi(x) = ; | pile = t) (ier =) a+ = | p(\ja —t|)(t — aj;-1) dt, Vi=1,...n—1. 


a Xi-1 


Then, the approximation equation of (1.22) has the form 
AUn (x) = S (wi(x)H (a, 2i)u(2i) + f(x), Vax € [a, b}. (1.23) 
i=0 
By choosing the collocation point 7;, we get the following algebraic system 


Mun(2j) = S “wi(aj)H (aj, 21)u(xs) + f(x;), g=0,1,...n. 
i=0 


We define forall x € [a,b] the error e,(x) by 


b n 
ata = / p( |e — t) (x, t)u(t) dt — Srun(2)H(@, 2i)u(a). (1.24) 
7 i=0 


Theorem 1.4.6. Let uy, be a solution of (1.23) and u the exact solution of (1.22). Then, un > u as 
n— oo 
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Proof. We have 


/\ 


b 
len(w)| [wlio | [ae t)u(e — tC, t)u(t dt, 


/\ 


b 
< / |p(\a — t])| a(e.t)u(e - = [a — t)A(a,t)u(a;_-1) + (€- vj-1) H(z, t)u(2:)| 
1 


+ i [(xs —t)(A(a,t) — H(a,2;-1)) u(ai_1) + (t — aj-1)(A(2,t) — (a, 2;)) u(zi)| dt, 


IA 


b b 
/ / t>( lor — #1)| [EE (2s) [4000 (ty A) + oo, o( EE, h) (a) ell. 


Then, 
enll: < | ma 2700) sh) ana (FF )(t) [ll laa 0-0). 


From (1.22) and (1.24) 


Iu — Unlloo S [lenI lx, 


and when n — oo, we get |len||x, where, ||en||x = sup |en(x)|. 
a<a<b 


Study of linear integral differential equations of the second kind: Analytical and numerical approach 


29 


Chapter 2 7 


Analytical and Numerical Study 
of Linear Fredholm 


Integro-difterential Equation 
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In this chapter, we begin the analysis of the integro-differential equation. We make some assumptions 
on the kernels to reformulate our equation to an equivalent system of integral equations. In section 
2.1, we study the existence and uniqueness of the solution. In section 2.2, we search for a numerical 
solution based on Nystré6m’s method, we follow some steps to show the convergence of the approximate 
solution among these steps the v-convergence of the Nystrém operator. In section 2.3, we apply two 
projection methods: Collocation and Kantorovich. In the convergence analysis, we show the convergent 
of the projection operator to identity operator. Afterwards, we have cried an iterative system of the 
three methods because the application of the last methods leads to a resolution of a very large algebraic 


system. 


2.1 Analytical study 


Let C1[a, b] be the Banach space of the continuous and differential functions with usual norm given by 


Va € [a,b], |lullerfa,s) = sup |v(x)|+ sup |v'(x)]. 
a<a<b a<a<b 
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We recall that, our linear Fredholm integro-differential equation is given as 


Va € [a,b], =f Ki (a, t)u wars | Kele,OW( dt + f(z). (2.1) 


Our problem in this section is to study the solution existence and uniqueness of linear Fredholm integro- 
differential equation (2.1) in the Banach space X. We need some conditions on the two kernels {Kp }i<p<2, 
that allow us to transform our integro-differential equation into a system of integral equations. Assuming 
that 


OK, 
H)| S*P(a,t) € 0%(a,02,R), p=1,2. 
Ox 
We obtain the possibility of defining the derivative u’ by 


Va € [a,b], Au! (x) = ee t)u(t) dt + es ? (a, t)u! (t) dt + f’(z). 


Now, we set up = u®-)) € Oa, b] for p = 1,2, then the equation (2.1) is equivalent to the following 
integral equations system 


Aur (2) - f Ki (a vt)ur(t yas Ko( (a ,t)ug(t ) dt + f(a iF 
Va € [a,¥ . (2.2) 
Au2(x y- [& aa (a, t)u,(t) dt + O(n, thualt) dt + f'(x). 


Its solution (u;,u2) must be found in the product space X, = X x X with the norm described by 
(1.1.11) where, X = C°[a, b]. 


Using the definitions and notations in 1.1.5, we define a family of linear operators {Tgp }i<q,p<2 from 
the system (4.8). For p= 1,2, we have 


b 

Va € [a,b], Vue X, Tipv(z) = / K,(a, t)u(t) dt, 

Va € [a,b],Vu EX, Topv(x) = as (x, t)u(t) dt, 
a Ox 


such that {Top }i<q,p<2 C BL(X). We represent Ar the block operator matrix in the following way 


Ar : Pe — Xe 


V ++ ArV(z) = (Tanto) + Ty2v2(x), To1v1(@) + Trata(2)), (2.3) 
Finally, the system (2.2) yields the next simplified version 


NU = APU +F. (2.4) 


The equation (2.4) has a unique solution if (AIz— A7)~' exists and is bounded as shows in the following 
theorem. 


Theorem 2.1.1. /f |A| > |||Av||| then (2.4) has a unique solution U in Xo. 
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Proof. We write, 


1 
(Alp — Ar)~* = Aa — yar). 


Then, we use Neumann’s theorem 1.1.6 the fact Hll4rlll < 1, to show that (Alz — Ar)~! exists and 


7 1 
lle - Ar)~* III < 


IAI — IArlll 


In the rest of our work, we need the compactness of the block matrix Ar. We state next two theorems 
to show that the block matrix Ay is compact. 


Theorem 2.1.2. All operators (Tqp)i<q,p<2 are compact in X. 


Proof. We prove that the subsets Tj,(Bx) = {Tgp(up), Up € Bxqp are relatively compact for all 
q,p = 1,2, where By = {v € X,||vu||x <1}. Firstly, we have 


Tip(tp)| <(b—a) Janes, [Kp(,t)| dt [hépllcojes, p= 12 


Tapl)] <(— od max, | Feo at lepllesoa. P= 12 


This proves that the sets T,,(Bx) are bounded. Now let z,y € [a,b] and u, € X. Then, 


ITiptip(#) —Tiptip(y)| < (b= a) max [Kp(@,t) — Kp(y,6)| [lupllx, 


OK, OK, 
[Topup() — Topup(y)| < (b—a) max ]—*(a, t) 


t : 
a<t<b Ox Ox (y, ) ||up| |x 


OK. 
Since K;,,(.,.) and ae .) are uniformly continuous, the sets Tj,(Bx) are equicontinuous subsets of X. 


By Arzela-Ascoli 2 Tj,)(Bx) are compact. Finally, we get that the operators {Tgp }i<¢,»<2 are compact 
in X. 


Theorem 2.1.3. The operator matrix Ap is compact in X5. 


Proof. Let the bounded set D = D, x Dy C Xo, such that V(z1, 22) € D,AC > 0, |\zi||x + |lzellx < C. 
We need to prove that the set 


Ar(D) = 1{ Stolen. He Dy} 


q=1 sp=1 
is relatively compact. For q = 1,2, we have that the sets 
2 
{Tonle Zp © Dp} , 
p=1 qd 
2 
are relatively compact because the operators T,, are compact operators and )> Ty, is also compact. 


p=1 
Therefore, 


{ST}. % € Dy} 


p=1 
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are a compact sets for g = 1,2, then by Tikhonov theorem [79] the product 


T1{ Stl) ZpE Dy} 


q=1 ‘p=1 


is compact. On the other hand, 


seb ee alte) in € Dy} 


q=1 ‘p=1 


| 
——— 
| 1S) 
oo 
8 
ea, 
er 
& 
B=} 
a 
S 
ae) 
Q 


q q=1 sp=1 


which proving that Ap(D) is relatively compact. 


2.2 Nystrom Method 


We recall A,,, for all n > 1 the uniform discertization of interval [a,b] given in (1.15) and apply the 
Nystrém method. We get, 


AU nl ya Lx; U1 nl n Xj) + Sure x 504 )Ug, al®) + f(z iF 
i=0 i=0 
Va € [a, b], 
OK2 
Au2 al =e (x, Xi) U1, yn (a) + ye * (a, Xj) U2,n (Xi) + f'(2), 


where the weights {w,;}”_, verifying 
n 


sup) el = =W <+o. 


n21; 
We select the collocation points 2;, 7 = 0,1,...n, to obtain the following algebraic system for all 
Méan(%j) = Doi Ki (xj, 4)urn(ai) +) wiKo (xj, 21)U2,n(2i) + f(z3), 


i=0 i=0 
j=0,1...n 


“OK 
Mutants) = Doi Ge (Bp Bt nl )+ ae 2 (25,2))uan(2i) + f'(23) 
i=0 


Now, we present the sequences approximation operators {Tyqn}1i<p,q<2, Vn > 1 and Vz € [a, }] 


Yue X, Tignv(z) = Soke xi)v(a;), g=1,2, 
i=l 
nm OK 

VUE Xx, Toqn¥(2) = 2a, x;)u(2i), q= 1, 2, 


and the Nystrém block operator matrix Ar, : a oe 


Tiin Thien 
Ar, = : (2.5) 
To1 n T22.n 
Finally, the Nystr6ém equation has this form 
Ol — Az, )UX =F, (2.6) 
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where, U' = (u1,n,U2,n)’ is the Nystrém solution. 


Convergence of the Nystrom Method 


The Nystr6m operator is known to be among operators who do not converge in the norm sense. In 
order to prove the convergence of the approximate solution to exact solution, we are going to show that 
our approximate operator is v-convergence. 


We explain what we have provided in this section in order to make it easier to follow the steps of 
convergence study. We start by proving that the Nystrém block operator Ar, is v-convergent to Ar. 
Then, we present the theorem which show that (Alz — Ar,)~' exists and is bounded. Finally, the 
approximate solution convergence is given in the last theorem of this part. 


In the next theorem, we show the pointwise convergence of the Nystrém block approximation. 


Theorem 2.2.1. Let Ar, be a Nystrom block operator matrix (2.5) and Ar (2.3) be a block operator 
matriz. Then, 


VUE Xo, jim ||(Ar, — Ar)Ullx, = 0. 


Proof. First, by the convergence of quadrature rule, we get 
Va € [a, }], lim [Tpqtq(@) — Tpgntq(t)|=9, 9,p = 1,2. 


We start by proving that our approximation operators een eer re are equicontinous. For each 
Uq € Oa, bj, x,y € [a,b]? and for q = 1,2, we have 


[Tigntq(#) — Tign(y)| < W max |Kq(x,t) — Kq(y, tI, 


a<t<b 
aK aK 
= < q q 
|ZoqnUq(@) — Toq,n(y)| ors (x,t) Oy (y, t) 


By the uniform continuity of K,(.,.) and 0K,(.,.), we deduce that {TpqnUq}1<q,p<2 are equicontinuous. 
Then, using lemma 1.1.1, it follows that all {Tp¢ntg}i<p,q<2 are uniformly converge to {Tpqv}i<p,q<2- 
Thus 


Vug € X, Jim |I(Zoa.n —Thq)Uq\lx =9, p,q=1,2. 


Finally, 


2 
(Az, - Ar)U||z = feces Flin — Tpq)qllx; 
a 


which gives ||(Az,, — Ar)U||x, + 0, when n — oo. 


KK K 
Let define the numerical integration errors for the integrand K,(,.)Kq(.,y), om (a, a, y) and 
x x 
Vu,y € [a,b], g=1,2 and Vn>1 as 
b n 
Eran(isy) =f Kale, QKg(t,y) dt — Yi Kq(v,2) Kyles), 
” i=1 


OK, OK, “. OK, aK, 
Eam(e) = fF Ao Gta) a Lore a) Fe). 
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For the following theorem, we prove that the above errors are uniform convergence to 0. 


Theorem 2.2.2. Letn > 1. Then, the errors {Epgn}i<p.q<2 are uniformly convergent to 0. 


Proof. We start by proving the equicontinity of errors (2.7). From the definition of equicontinuty 


V (x1, 22) € [a,b], Ve> 0, a0 0, V(yi, Y2) € [a, b}°, ||(21, 2) _ (yi, Y2)||R2 < 6, 
we need to prove that 
\Epgn(21, y1) — Engn (La, Y2)| < E, G,p= 1, 2. 


We have, for g=1,2 and p= 1 


lErg.al@as y1) = Eig n(2, y2)| < |Erqmn (2x, Yi) a f1ge (01, Ya) + |E1qn(©1, Y2) _ E1qn(2, Y2)|; 


with, 


[E1gn(21, 91) — Ergn(21,y2)| < (b-a+W) max | Kq(«1,t) — K,(x2,t)|, 


|E1gn(21, Y2) — Ergn(r2,y2)| < (b-a+W) max | k(t, yi) — Ka(t, y2)|- 


In the same way, for p = 2 we have 


exerts y1) = Exq n (2, y2)| < |Ezqn(21, ¥1) a E2q,n(1, Y2)| + |E2q,n(£1, Y2) ol E2q,n(2, Y2)|; 


with 
OK OK 
|E2q,n(t1, ¥1) — Ergn(t1,y2)| < (b-a+W) es aa (x1, t) ae (x2, t) 
OK OK, 
|E2q,n(#1, yo) — Ergn(t2,y2)| < (b-a+W) eee an (t, y1) a (t, y2)|, 


OK, 

The Hein theorem’s [48] proved that Ky and oe are uniform continuous in the compact set [a, b]?. 
xv 

Then, 


ign (Gist) = ig nltesgs)| =< 22, 


|Eoq,n(@1, 41) — Erqn(@2,y2)| < 2e. 


By Lemma 1.1.1, we get that {Ep9n}1<p,q<2 are uniformly convergent to 0. 


Theorem 2.2.3. Let Ar and Ar, be the block operator matrix present by (2.3) and (2.5) respectively. 
Then, 


(Ar —Ar,)Ar, ||| = 0, and (Ap — Ar,) Ar, lll = 0. (2.7) 


lim ||| lim ||| 
n— +00 n— +00 
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Proof. For n large enough, we have 


Ti - Thin Tie -— Tian Ty Tig 
(Ar —Ar,)Ar = 
To} — Tain To2 — T22,.n To, Toe 


Ty _ Thin Ti2 _ Ti2,.n T, sn T12,n 
(Ar—Ar,)Ar, = 


To — Tain To2 — T22,n Tein Toon 
Using Hadamart product (product element by element), we get 


Ti=—Tinta Cs—Dien)Tie 
(Arp — Ar, )Ar = 
(T21 — Toin)To1  (T22 — Te2,n)T22 


(Ty =TigTaa (Te Tien) Tan 
(Ar —Ar,)Ar, = 
(To1 — Toin)Tain (Tee — T22,n)Tr2.n 


By the norm of block operator matrix, we obtain 


(Ar — Az, )Ar||| = Pere Yee pq ~ Lpqn)Z pall, 
S182 5 
II(Ar — Ar,)4r, [I] = a, De pa — Tran) Tpaanll 


where V 2(@-)) € ©°[a, b] and q = 1,2 


b 
(Tiq = Tign)Tig2—» (a) = | Ergn(z,y) z'-U(y) dy, 
b 
(Taq = Toqn)Toqg2"4») (a) = / Erqn(Z, ¥) z'a-U(y) dy, 
(Tq — Tign)Tign2 (a) = yw Erq,n(#, 24) 2) (a4), 
(Tag Then) Pant) = So Egan) 2"), 


I|(Tiq ~ Tign)Ti4l| 


IA 


(b — a) max JEign(a, y) F 


(Pia ~ Than)Pian|| SW max lEign(a,¥)|, 


a<x,y 


| (Taq = Toq,n)T2q\ | 


IA 


(6—a) max |Eqn(2,y)I, 
anys 


|(Loq — Taqn)Toqnl| < W = |E2q,n(z,y)|- 
a<ax,y<b 


Finally, by theorem 2.2.2, it follows that max. |Epq.n(, y)| 7 0asn- oo. 
asna,ys 
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2.2. Nystrom Method 


Theorem 2.2.4. Let Ar, be a Nystrém block operator. Then, (AIz — Ar,)~' exists and bounded. 


Proof. For large n, we have 


Ig + (Alg _ Ar)~' Ar, (Alg _ Ar, ) — (Ala = Ar)! [Ar — Apt Ar, 5 


n 


(Aly _ Ar)! or Ar, ) MAr Ar ) 


+ (Ag- Ar, Ar] ; 


= (Alp— Ar)! are — Ar)+ (Ar - Ar,)Ar}, 


Ng + (Alp — Ar) ' (Ar — Ar, )Ar- 


Then, 

1 1 

X In + (Aly - Ar) ‘Ar, (Aly - Ar, ) =Ig+ x Ale - Ar) '(Ar - Ar, )Ar. 
In the other side, we have 


1 


1 
|< (Al, — Ar)" (Ar — Ar, )Arll| < 
r |A| 


II(AZ2 — Ar)" I(Ar — Az, )Arll]. 


But, in the previous theorem |||(Ar — Ar,)Ar||| 4 0 when n > oo, we get 
1 


|< (Alp — Ar)" (Ar - Ar, )Arl| <1, 
Xr 


1 -1 
which proves that by Neumann’s theorem | [2+ rl (Alp — Ar) (Ar — Ar, )4r) exists and is bounded. 


Then, the inverse of (AIz — Ar,) has the following form: 


-1 
(Alp — Ar, )7' = ( + (Alp — Ar)7'(Ar — Ar,)Ar € + (Alg — An)" Ara , 


and 


-1 


WAL2- Ar) II] S NV+ Qh - Ary" Ar ~ Ar,) Ar) III Zz + (Ala — Ar)~* Ar, III, 


Z II|Z2 + Ale = Ar)~" Ar, |II 
~ || [Ila A) Ag = Ar, Ag |||" 


Finally, we give a majoration of the error in the next theorem 


Theorem 2.2.5. Let UN = (uin,U2n)! be the Nystrém solution of (2.6) and U = (ui, u2)* be the exact 
solution of (2.4). Then 


Jim ||Un — Ullx, = 0. 
Proof. From to equation (2.6) and (2.4), we have 
MUn —U) = Ar,UN — Ap = Ar, (U — UN) 4+ (Ar, — Ar)U. 
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Then, 
eal -U= (Aly _ Ar,)7* (ar, _ Ano | : 


We get 


On’ — Uz, S MWOALe — Arn) "II MI(Ar, — Ar)U IIx, 


and when n — oo, we get the result. 


2.3. Projection Approximation 


In this section, we give the basic concept of projection method see [5, 2, 48]. In order to illustrate 
the general framework of Collocation and Kantorovich methods. First, we recall A, (1.15) the uniform 
discretization of the interval [a,b]. We denote by Xan the sequence of finite dimension subspace of the 
Banach space X2. The principle of both methods is based on the application of projection operators 
{Pr}n>1 defined by 


PR, : Xa — Non 


V +>P,V(2)= (Finer(v), Ponta(e)), (2.8) 


such that {Pynv}n>1 for g = 1,2 are sequences of projection interpolation. This can be written in the 
following linear combination of piecewise linear functions. 


Yn >1,Va € [a,b], Pynvg(x) = S“aiei(x), 
i=0 


where {a;}"_) are unknown coefficient to be determined and {e;}"_, are the hat functions. These are 
given by (??), such that {P,}n>1 satisfies the following condition interpolation 


PonVq(ti) = Ug(ai), t=0,1...n. 


_ We start by showing that our projection is pointwise convergent to the identity operand in the space 
Xo. 


In the following theorem, we show the convergence of P,,U to U in Ca and then, the convergent of 
Ar to P,Ar. 


Theorem 2.3.1. Let U = (ui, uz)’ € X2 and {Py}n>1 be a projection operators sequence given by (2.8). 
Then, 


[|(Z2 — Pa)U| |x, < t00,1(U, h), 


and when h + 0, we get ||(2 — Pn)U||x, + 9. 


Proof. For n > 1, it is clear that 
I|Z2 — Pr)U |x, = Joes eae) — Pi ntii(@)| + max, |ua(2) — Pente(z)|. 


We have, for all (ui, uz) € X2 
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Jui (x) — Py nta(x)| = |u(x) - do (i) ex(@)| 


= S “ei(2) [ua(2) — ui (z%)| 


1=0 
= Diet )| Jus (x) — ur (as) 
< Koo,0(U1, hi); (2.9) 


such that, | >> e;(x)| = 1. In the same way, we get 
i=0 
|u2(x) — Penu2(z)| <  Koo,o(ua,h). (2.10) 
Substituting (2.9) and (2.10) in (2.9), we obtain 


|(Z2 — Pr)U|| x, S Koo,o(t1, h) + Koo,0(Ua, h). 


Lemma 2.3.1. Let Ar be bounded and compact block operator matrix and {Pr}n>1 converges pointwise 
to the identity operator. Then, 


‘lim |[|(Z2— Pa) rll = 0. 


Proof. Since Ar is compact 


M = {ArU,||U|lz, < 1} 


is relatively compact sets in the Banach space X and by the Banach Steinhaus theorem {Pa}nen con- 
verges uniformly to Jj in M, i.e 


lim |||(Z2 — Pn) Ar||| = Jim sup ||(I2 — P,)ArU |x, = lim sup ||(Jo — Pr) 2 | lye 
noo ONT |g NCO FEM 


In the numerical treatment, we need to verify the existence of (AIz — P,,Ar)~} to demonstrate the 
convergence of the approximate solution u, to our solution. We introduce the following theorem to 
ensure this existence. 


Theorem 2.3.2. Let P, Ar be a projection approximation of block matrix operator and Ar given by 
(2.3). Then, (AIg — P, Ar)~! exists and bounded. 


Proof. We have 


(Alp — Pn Ar) = (Alp — Ar) [Io — (Al2 — Ar)7 (Ar — Pr Ar)). 


Using lemma2.3, we get 
Jim |||(f2 — Pa) Arrll| = 0. 
This shows 
[| Alan — Ar)~"(Ar — PrAr) ||| < |[]QL2 — Ar)“"III [I|Z2 — Pa) Arlll < en, 
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which gives 


I Ala — Ar)? (Ar — PrAr)||| <1, 
so that 
Zi 
(1 = (Alg = Ar) \(Ar = PaAr)) ’ 


exists and 


| (te- (A - Ar) P,)An)) 


Finally, (AIg — P, Ar)~' exists and |||(AIz — P, Ar)~||| < C, where 


WAL — Ar) "Il 
7 1l—en , 


C 


Now, we present two projection methods. The first is Collocation method which is introduced by 
Atkinson [5] and Nair [48] for searching the sequence US € Xo C X2 solution of 


MUS = P, ApUS + PyF, (00) 


whereas the second one is for finding the sequence ue OnE Xan S XG solution of 


AUS = PA ee (2.12) 
It’s the Kantorovich method. 


2.3.1 Collocation Method 


The main idea of collocation method is to find a sequence UC solution of 
AUS = Pr ATUL + PrF, (2.13) 


where P,F = (Praf,Paf’)', and Un = (uf, uf,)* is Collocation solution are given by 
Urn(e) = de avei(x), 


uf,(2) = Y Bres(a). 


where {a;}o<i<n and {G;}o<i<n are coefficients to be determined. 


The next theorem shows that the approximate solution US converges to the exact solution U. 


Theorem 2.3.3. Let U = (ui, ug)! the solution of (2.4) and US is the Collocation approximation of 
(2.13). Then, 


IU — Un Ilx2 SC] oof, h) + [Iz — Pr)Arlll]- 
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Proof. We have, for n large enough 


U—US = (Ab — P, Ar) E —P,F+4+(Ar- PAnyU). 


Then, 
IU —UEllg, < lb — Par) I — P,)Fllz, +|\(Ar — PrAr Ulla) 
< oli ~ P,)Fllg, + Ill — Px) rll: 


< Coal h) + [IIo — P,)Anlll): 
Also, we have 


I|U — OF Nx = [ttn = lle fa,s): 


Then, when n — 00, we obtain uf, = uf > u in C1[a, b]. 


System Approximation 
The system (2.13) can be written as 
dul» (2) = Pr Tigtty 4 (®) + PF igs (a) + Pi f(x), 
Va € [a, Bl, (2.14) 
Aus, (2) = PnT21U{ n(x) + PrT22U§ n(2) +P, f' (2). 


These are equivalent to 


Sayes(e) = So rufa) + Tigu,, (4) + i) a: 
J= I= 
> Byes(a) = > [tavul (es) + Thouf,, (xj) + F(a) e;(2). 
I= = 


This leads to the resolution of following block matrix system form: 


aj = An G, ta + 2d Aral3; i) Bi +f 


i=0 


py = 2 Aoi (J, ta + 2 Arald; i)Bi + fj, 


where f; = f(x;), fj = f’(x;) and for p= 1,2 
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1 Liq 
Aip(j,t) = ; | K,(«;,t)(t — vi-1) a+ | K,(2;,t)(tiz1 — #) dt, O< jg <n,1<i<n-l, 
Tj-1 Ly 
1: of 
Ap0) = Gf Kyle ler) dt, OSI <n, 
a) 
a i . 
Aip(j,n) _ a Kglty, tit— ty) dt, O<j<n, 
Ln—-1 
_ 1 (” OK Tl OK ; . 
Aap(5,t) = aA pacleavtN(t— aia) ae f Hey (tir (wits — t) dt,0<j<nil<i<n-l, 
1 f® OK . 
Amp(0) = 5 f FRester-t dt, O<F<n, 
xo 
1 f*" OK 
Arp (j.n) = =f Hy (tar tlt — @n—1) dt, O0<j<nl<i<n-1. 
gu OF 


2.3.2 Kantorovich Method 


In the collocation method, the convergence of the approximate sequence is guaranteed by the following 
assumption P,, f > f, when n > oo. We can cancel this hypothesis, if we assume that (f, f’) € Ran(Ar). 
For this raison, we apply on both sides of the system (2.4) the block operator matrix Ar to get 


AAU = ApArU + Arf, (2.15) 
We denote 
Ty The u 
V=A7U = ; 
To, Too u’ 
and (2.15) is equivalent to 
AV = ArV+ ArF. (2.16) 


If V = (v1, v2)! is the solution of (2.16), then the system (2.4) can be rewritten as 
Kan 1 


We compute an approximate solution in the following form: 


1 
or = = (ve + F), (2.17) 


such that V,C = (vf,,v§/,)! is the Collocation solution of 


MW,C = P, ArV," + P,ArF, 
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represented in the following way 


ve, (e) = dare), 
Va € [a,b], 7 
V5 n(£) = > Biei(2). 


Therefore, we can define Kantorovich’s approximation as 
BB oat = PAp a FP, (2.18) 


with Ue" = (ufiAn, uger)*. We notice that from the two equations (2.18) and (2.17) 


an 1 
|r" —Ullx, < pylYar — Vall. 


Now, we state a theorem which shows the Kantorovich approximation (2.17) convergence to the exact 
solution. 


Theorem 2.3.4. Let U = (ui,u2)* be solution of (2.4) and US" = (uke, usin")? is the approximate 
solution given by (2.18). Then, 


| — Ue" lle, SC |I|a — Pr) Arlll: 


Proof. For n large enough, we have 
U —Ue = (il — P,Ar) (F - (Ah -— Ar) F 


= (Alp — P, Ap)~* E — (Al, — P, Ar)(Alp - Ar) 1F] 


= (Ah — P,Ar) lian Ar)U — (Al PyAn)O| 
= (Ah — P,Ar)" [Are - PyArU| 


Then 


| — US" Iz SC ||A2 — Pa)ATU lx, 


< C |||(e — Pr) Ar||. 


This theorem shows the convergence of Kantorovich approximation solution to the exact solution 
when n — co in C'[a,b] because 


| — Oe lhe, = Ne = wr lot fab): 
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System Approximation 
Applying the Kantorovich method, we get the following block matrix system 
n 


aj = DA, OG + Aral i)B; + Gi(j), 


i=0 


By = S7 An (j,i)or + 3 Anal, i) + Gas), 


1=0 i=0 


such that for 7 = 0,1,...n 


for p= 1,2 and0 <j <n, we have 


1 Hie itd 
Aip(j,4) = a Kp(zj,t)(t — 24-1) a+ | Ky (aj, #)(te41 —t) dt, 1<i<n-1, 
; 1 
Aip(i.0) = 5 f Kylej.t)(0r 2) a 
B20) 
| 1 ft 
Aip(j,n) = ; Kp(xj,t)(t — tn-1) dt, 
_ 1 [* OK "+1 OK. . 
Agy(j,4) = ae. on a ate ff Fy (tart) (wisa — t) dt, 1<i<n—1, 
: 1 ("OK 
Anp(i0) = 5 f  Z¥lej t)(er-8) a, 
xo 
; 1 f*" OK 
Anglin) = 5 f° F(a; t\(t- ena) at 


b b 
[Rules .40) a+ | Kola; f'(@ dt, §=0,1,..., 


i K 
- [{& w ——*(2;,t) riya f —(x,;,t)f/(t) dt, 7 =0,1,...n. 


2.4 Projection block iterative scheme 


From the previous results, we understand that if we have n too large, we get a good precision on the 
solution and a good error convergence to zero. But, it poses problem because the size of the block matrix 
P, Ar becomes large. Therefore, the conditioning of the matrix can be affected. An iterative scheme for 
solving the two projection schemes presented by (2.11) and (2.12) is necessary. In [11], authors proposed 
a numerical method for solving 2 x 2 block systems and may be applied to solve this system especially if 
the coefficient matrices are ill-conditioned. In this section, we use the block iterative scheme in order to 


avoid calculating the inverse of the whole matrix (AIz — P, Ar)~! 
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So, for p = 1,2, |A| > ||Tpp||, then (AI —T)»)~! exists and bounded. We can present the matrix (AJz— Ar) 
in the following form 


Ala — Ar 


AI -— Ty O I O O (AL = Ti) 'Ti2 
O AT — To oO ! (AL _ To2)1T21 O 


The main idea of this scheme comes from (Aly — Ar) introduced by (2.19) and used to construct Collo- 
cation block iterative scheme 


AI — PrT 1 O Y Paf 
O AL —_ PrTy2 Y2 Bf’ 
where 
Y1 I =F =P, Fi)" Pate ne 
Y2 —(AI — Py To2)~* Py T21 I ulG 


Therefore, we write the Collocation scheme in the following form for all n > 1 


yr = (Alg — PaTii1)7*Prf, 


yo = (Alo — PrT22)7'Prf’. 


CO 
Un =Y1, 

CO 
Ugn = Y2; 


ce = (Ala — Py T11)~1 PaTizuy 3" T Yi; Vn,k = 1, 


k as k- 
tig, = (Alp — PyT22) 1p, Tour thy Wn,k> 1. 
With the same principle given above, we build the Kantrovich iterative scheme: 


zy = (Alo — PaTi1) 71 (PaTuf + PrTiaf’), 


29 = (Alp — PrT22)7!(PrToif + PrToof’). 
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2.4. Projection block iterative scheme 


Kan,O __ 
1,n = 1; 
Kan,O  __ 
2,n = £2; 


ue = (Aly _ Pt aa vT 41; Vv n, k 2 1, 


1,n 


Kan,k _ kan,k—1 
linn == (Alg— P,Ta2)* Patou, +z, Vn,k>1. 


In the next theorem, we prove that the iterative scheme of Kantorovich and Collocation methods 
converges to the exact solution. 


Theorem 2.4.1. Let UF = Che ee V k,n > 1 the iterative solution of projection method and 
U = (ui, ua)’ solution of (2.4). Then, 


. k 


Proof. We have 


\|U — Urllx, < MIU —Uallx, + llUn — Ulex, 


such that 
O (AI — PnT11)7*PrTi2 Un — Urn 
U,—-uk = 
(AI — P,T22)7pnToi O U2,.n — i 
By recurrence, we get 
O (AI — PaTi1)PpTia\ * (urn — Uo np 
Uy —U® = 
(AL = PyT22)~!PnTa1 O U2,.n — 0 
We denote Mp, r,, the block matrix given as 
O (AL = Ptah “PTs 
MP. Ton = 
(AI — P,T22)~*P,To1 O 
Then, 
Un ~ US = MB.r,,(Un- U2), (2.19) 
where 
Uy — UL = Mryn,, (UU +0). (2.20) 
Substituting (2.20) in (2.19), we obtain 
[Un — UEllg, € I Mese lll [Ikon — Ullz, + UCI. (2.21) 
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Then, by the inequality (2.21) 
7 —Unllz, < UO -Uallz, + Un — Valle, 


< (1 [lates rap! ) IO — Gall, + Mest ICL gy 
Now, we will prove that |||Mp,7,,||| <1. By usingthe norm given in (1.1),we get 
HM, qll = max (AZ ~ Py Taa)~“Py Fal (AE ~ PaTin) Pal). 
But, 
(AI — P,T22)~*PrToi|| < ||(AL — PrT22)7"|| || PrTaill- 
Also, we have 


AZ — PrT22)* || << ||( — (AL — PaT22)7*(T22 — PrT22)|| |AT — T22) "II, 


. (AT — Ta)“ 
= T=] Ts) — Pa Daal 


Theorem 2.1.1 shows that Jim |i — P,,)To2|| = 0, we get 
||(AZ = PnT22)~*|| < ||(AT — T22)7"|I- 
Then, 
||(AZ — Pn T22)~* Px Tail] < |Z — T22)7*| ||T2il| < 1. 
In the same way, we get 


(AZ = PaTi1)~* PaTial| < ||AT - Tir)" || ||Zial| < 1. 


So, |||Mp,7,,||| < 1 which proves that 


. k 
Yn =A, lim || — UR l|x, < IO — Tall x 


We prove in theorem 2.3.3 and 2.3.4, 


Jim || — Ualls,. 
which gives 


2 2 _ 7k = 
Jim (im IU) = 
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2.5 Numerical Tests 


In order to show the effectiveness of the suggested methods, we apply these approaches to the following 
numerical examples. We use the following estimate to calculate the error between the exact solution and 
the approximate solution at points x; 


err = max |u(x;) — un (as)| + max lu'(as) — up (a4)L, 


b-a 


where the points x; = a+ th, forh= andi=0,1,...n. 


nr 


2.5.1 Test O1 


We consider the following integro-differential equation 


Vx € [0,1], Au(x) = | * 4) ae | sie f(z). (2.22) 


ev + ef 


We choose \ = 2 and the exact solution u(x) = e® to get: 


1 arctan s 
f(x) = log(e® + 1) — log(e! + e”) + 2e* arctan ( c ) ( aa) 


Table 2.1: The error between the exact and approximation solution of equation (2.22) 


n Nystr6m Collocation Kantorovich 
10 = 1.6368e-04  6.6852e-04 1.0445e-05 
100 = 1.6355e-06 —_7.1682e-06 1.1061e-07 


-4 
7 x10 . 
—— Error of Kantorovich metod 
—+*— Error of Collocation method 
6f --£+- Error of Nystrom method 
5 bese 
AL J 
37 J 
27 4}----8----4}----g | 
. A i Gel ini ia Cee 
Senge 
i J 


Figure 2.1: Errors of Nystrém, Galerkin and Kantorovich Solutions of (2.22) with n=10. 
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2.5.2 Test 02 


We consider the following integro-differential equation 


; ate lt eee 
@eh+e2ei" [e@ + A)! (t)dt + f(a). (2.23) 


Va € [—2, 2], \u(x) | 


We choose \ = 4 and the exact solution u(x) = x”, to obtain 
2_ 1 2 2 
f(a) = 4a* — 3 (arctan(e + 16) — arctan(zx )). 


Table 2.2: The error between the exact and approximate solution of equation (?7) 


n Nystr6m Collocation Kantorovich 
10 0.2315 1.4478e-03 1.6786e-05 
100 =: 0.0034 1.7149e-05 1.8614e-07 


0.25 T if T T T T 
—— Error of Kantorovich metod 
—+*— Error of Collocation method 
a --4-- Error of Nystrom method 
\ 
0.2 L ‘i \ Fi 3 ma 
Fe x : ‘ 
¥ \ 3 
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/ \ 4 tk 
/ \ ‘ 
0.15 ¢ / \ ' \ 1 
/ \ 
7 \ | x 
\ ry 
a \ 4 \ 
/ Pe Ed \ 
/ a a \ 
0.1 + / “a \ 7 
F \ 
/ \ 
/ \ 
/ q 
/ \ 
/ L 
/ \ 
| r % il) 
0.05 a a 
a” ‘a 
03—_—-—_ + + % ¥ $e H 
-2 -1.5 -1 -0.5 0 0.5 1 1:5 2 


Figure 2.2: Errors of Nystrém, Galerkin and Kantorovich Solutions of (2.23) with n=10. 
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2.5.3 Test 03 


We give third equation in the following form 


10 u 10 au! 
Var € [1,10], Au(a) =| Ss a+ f an dt + f(a). (2.24) 


We choose \ = 20 and the exact solution u(x) = log(x), we get 


f(x) = 20log(x) — 2\/(a + log(10)) + 2.x — arctan(10log(10) — 10 + x) + arctan(x — 1). 


Table 2.3: The error between the exact and approximation solution of equation (2.22) 


n Nystr6m Collocation Kantorovich 

50 =.2.3856e-04 ~—-6.7363e-05 8.6247e-06 
200 =1.5078e-05  4.2353e-06 5.3879e-07 
500 2.4140e-06 —6.7791e-07 8.6205e-08 


-4 
25 * : 0 T T T T T T T T 
—— Error of Kantorovich metod 
—+*— Error of Collocation method 
--4-- Error of Nystrom method 


Figure 2.3: Error of Nystrém, Galerkin and Kantorovich Solutions of (2.24) with n=50. 
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Chapter 3 7 


Analytical and Numerical Study 
of Linear Fredholm 
Integro-differential Equation 


with Weakly Singular Kernels 
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Recall the aim of this work is to make a complete and thorough study of the integro-differential Fred- 
holm equation. Therefore, we would like to focus on another form whose equation takes. Thus, our 
objective in this part is to study the possibility that both kernels of our equation are weakly singular. 
For this reason, we start in the section 3.1 with the construction of a sufficient condition that gives the 
existence and the uniqueness of the solution. Subsequently, we propose two different numerical methods 
which are: Product integration method in section 3.3 and Cubic collocation b-spline method in section 
3.4. In order to find the best numerical solution and reduce the calculation time. We present some theo- 
rems to ensure the convergence of the approximate solution. In section 3.5, we give numerical examples 
to see the error behaviour of two approximate methods. 


In this chapter, we study the following linear Fredhom integro-differential equation in the Banach space 
C*[a, 0] 


b b 
Va € [a,b], Au(a) = / pi(\x — t|) Ki (a, t)u(t) dt +f po(|x — t|)Ko(ax, t)u'(t) dt + f(x), (3.1) 


where \ #0 is a real or complex parameter, the kernels K; verify the next assumptions for i = 1,2 
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OK; 


Fa (x,t) € C%((a, b]?, R), 


M;>0, max (Ixi(e.0), 
a<az,t<b 


OK; 

a (2,0)| | < Mi, 
me (est}|) < 
and the singularity part p; verifies the following assumptions for 7 = 1,2 


PE W140, b—a). 


cf / — 
lim, |pi(s)] = +00. 
For 7 = 1,2, p; are chosen to take one of the following forms 


1. p(x) = |x —t|**1, a eE)0,1[, 


2. p(x) = |x — t/log(|a — t|) — |a — tl]. 


Under the assumptions (H1) and (H2), we can give the derivative u’ implicitly as 


OK, 
Ox 


Va € [a,b], Au’ (x) = / pi(|z — t|) (x, t)u(t) dt +f sign(«x — t)p} (|e — t|) Ki(a, t)u(t) dt 


+ / pa(|a— t)) 2 (x, t)u'(t) a+ | sign(ax — t)p>(|a — t|)Ko(a, t)u'(t) dt+ f’(x), (3.2) 


where 


. 1, if «>t, 
sign(z — t) = { —1 else. 


We cannot move on to the search for a numerical solution without proving the existence and uniqueness 
of the exact solution. This is done in the next section. 


3.1 Analytical Study 


The aim of this part, is to use the hypotheses mentioned above to construct a sufficient condition 
that ensures the solution existence and uniqueness of the equation (3.1). So, let us define T the linear 
integro-differential operator by 


T: C [a,b] —> C"[a,d 


b b 
u > Tu(2) -| pi(|x — t|) Ki (a, t)u(t) asf po(|a — t|)Ko(a, t)u'(t) dt. 
Thus, the equation (3.1) has an equivalent form, which is given by 
(Al —T)u= f, (3.3) 


where J is the identity operator of the Banach space C"{a, b}. 
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To prove that the equation (3.1) has a unique solution, we need to prove that the inverse of AJ — T 
exists and is bounded. 


Theorem 3.1.1. The linear Fredholm integro-differentail equation (3.1) has a unique solution u € C*[a, b} 
if 


Jal > 2] Mallellvsop-n) + Malealss (oso) 


Proof. We have, 


b b 
[Tu(z)| < / Ipi(|x — t])| |i (a, t)| |u(t)| a+ [ |po(|x — t|)| |Ao(a, t)| |u’(#)| dt, 

<  My||pil|r240,e—a |!ul|cofa,e) + Me||pellr2jo,e—a) IIe’ Icofa,0}- 

and 

/ Ok, : / 
|(Tu)'(x)| < Pd |x — t])| Fe tt) |u(t)| dt+ [| |pi (|x — t|)| |Ki(2, t)| |u(t)| dé 
OK2 : ! / 

4 ak |a — 2))| lpg t)| |u'(t)| dt+ | |po(x — t|)| |Ko(a, t)||u’()| dt, 
< My||p1||w2.2J0,e—¢) [|u| |eofa,o) + Me||p2||w12[0,0—) |!u"||cofa,t): 

Then, 


Tullotta.s) < 2] Maloalhvs ano + Mallpallvso-a) |lullera 


which gives , 
E21 < 2] Malipallsagoa-o + Mallaliwsayoo-c | 
Using the fact that 
|A| > 2) Af leulhv oa + Ma||pa\|w21[0,0—a] IF 
we get 
ITI < |Al- 


Therefore, according to theorem 1.2.11, we have sp(T’) C B(0,||T||), where B is the ball center 0 and 
with radius ||T7'||. Therefore, 


fr € C*(R*), |A] > 2[Mil|pil|w22[0,0-a) + Me||p2Ilw21 [00-4] i} C re(T). 


By Neumann’s theorem [5, 2, 89], (AJ —- T)~! and 


1 


|AL-T) "|< 
|A} — 2[Mi||p1||w21[0,o—a] + Ma||p2||w22[0,0—a)| 


This proves that the equation (3.1) has a unique solution. 
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3.2. From the equation to the system 


In the numerical treatment, we start by the construction an approximate method based on product 
integration. 


3.2. From the equation to the system 


In the analysis of convergence of the product integration method, we need to prove some proprieties 
on the linear operator T’.. So, we transform our operator T to a matrix of block linear operators. We put 
X = Oa, b] and X, = X x X. 

Let defined the projection operator P and the injection J by 


P: X_ — Ca, d 


(v1, V2) r— P(v1,v2) = Pv, 
J: Y ~ +X 
v + J(v) = (v1, v2), 


We denote the function u'?~1) by u, for p= 1,2. Then, 
T = PArJ, 
with Av is a block operator matrix defined on X> into itself by 
Ty The 
Ar = ’ 
Ta The 


such that, for r = 1,2 


Va € [a,b], Vu € X, Ty,v(x) 


b 
ii p,(\x — t|)K,(a, t)u(t) dt, 


b 


b 
[ volie— Gotu ats [ sign(x — pl ir =e) Kee, to(t at 


Va € [a,b], Vu € X, To,v(x) 


and the norm of Ar is defined by 


2 
| Ar ||| = max SollZirll- 
srs 


Finally, the equation (3.3) has an equivalent form is given by 


NU = ApU+F, (3.4) 


where U = (uy, ug) et F = (fi, fo). 


3.3 Product Integration Method 


Let n > 1. We recall that A, is a discretization of the interval [a,b], which is given by (1.15). By using 
the product trapezoidal rule 1, we get for all t € [xj;_1, 2], ur € X and r = 1,2, 


rei Oe = (03 (ee aides (i ase ae, 
SS (x, duel] = (03 - 1) (a, De 4 Ue) “(i - 25-1) Se" (w,j1 Wel 3) 
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Then, we can introduce the approximate operators for all v € X, x € [a,b] and r = 1,2 as 


Tipnv(t) = S “wr, j(2) K(x, 2;)u(23), 
j=0 
n aK, n - 
Tornv(t) = vrs) Z-(@, x;)v(aj) +) @,,j(x)sign(« — xj) K(x, x;)v(x;), 
7=0 j=0 
with, 
1 /™ 
Uro = a (a1 — t)p,(|a — t]) dt, 
Xo 
1 f* 
Urn = 7 (t — @p_-1)pr(|x — t|) dé, 
Ln-1 
1 tits 1 
wry =f (ep Opelle tl) ate =f” (b= 2y-1)pe( le =A) dt, 1 m1, 
x5 Uj-1 
and 


a Lp 
Dro = Zp er Oeelle— th at 
zo 


= Le 
inn =f (t= anaa)ph(le= tl) dt 
Bn-1 


no 


Tj41 1 f* : 
wing = Gf Ophea) ate f(t a5-a)po(le—t)) at, 1S j<n=1. 
x Lj-1 


j Paes 


So, we get an approximation system of the following form 


Arn) = YE wrg(0)Ba (3) n (85) + So weg) Kole, 25) an), 


Auan(z) = ¥ wns) GA Lj )Uin(L;j) + DS Wei (x)sign(a — v5) Ky (x, £;)Uin (a,j) + fi(x) 
j= 


ss 
iM 


n fa) n 
+ Dw rg 2) Ew, 2j)e u(y) + Brg a) sign ~ 25) (2,23) 02(05) + fol2), 
j=0 j=0 


By using the collocation points x = x; for i= 0,1,...n, we get the algebraic system: 


Aurn(zi) = D0 weg (ai) K1 (21, 23) Ua n(x) + D0 wr,j (ei) Ko (ai, £3) Uan(aj) + f(z), 
j=0 j=0 


n OK n 7 ; 
Nu2n(@i) = De wr (wi) (wi ty) n(@5) a 2 tor, (ts) sign (ae — ©5) Ky (xi, 27) Un (a5) 
j= j= 


n OK n 7 ; 
+2 wg i) Be (is 5) U1 n(5) + di tir, 5 (wa) sign(a = By) Ko (as, Gs) (ae) + Pole), 
j= j= 


which is equivalent to the next matrix equation 
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Mp Ante eo 
where U, = (tin, U2,n) and 
Ti1,n Ti2,n 
Ar, = 
Ty n Too n 


Convergence Analysis 


The convergence analysis of the product integration method is based on essential step. So, we start by 
presenting a theorem and a lemma. We follow all the steps that are contained in it. 


To show the following convergence 


||U — Unllz, 4 0, when n > oo, 


we will demonstrate that Ar, verifies the following steps in order to apply theorems 1.1.8 and 1.1.7. 


Let Ar be a block operator matrix defined on the Banach space X9 into itself and {Ar, }n>1 be a 
sequence of an approximation operator Avy. Then, 


(i) The operators Ar and Ar, are linear (clear), 


(ii) The block matrix operator Ar,, for n > 1, verifies this pointwise convergence 


We X2, lim ||ArU — Ar,U|]z, = 0. 


(iii) The set {A7,,,n > 1} is collectively compact, which means that the set 
5 = { ar,0,n > 1,|WUils, <1}, 


has a compact closure in X9. 
We present the following theorem and on which we are going to use the above mentioned steps. 


Theorem 3.3.1. Let {Tipn}i<ir<2 be an approximation operator of {Tir}i<ir<2 . Then, the sets 


{Tien Ola, Hh ltallodag Saws i} 


1<1,r<2 
are pointwise bounded and equicontinuous. 


Proof. Let uy € X and x € [a,b]. For r= 1,2, we have 


ITirntr(2)| < M,||Pr||71[0,o—a]||Ur||x, 
[Torn Ur(2x)| < M,||p-||12[0,0—a] ||tr| |x, 
which gives, 
|Zinnl| << Mr||pr||i+{0,r—a); 


[Zirmll << Mrlpylle+(o,e-a)- 
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This implies the uniform boundedness of T),,. From the Banach-Steinaus theorem 1.1.3, T),,, are point- 
wise bounded. Let prove now the equicontinuous. For every ¢ > 0 there exist, 6 > 0 such that for 
x,y € [a,b], |x —t| < 6, we have 


is eld Ta teh: < / bpe(|e — tl) [Kr(w, tu (#)]n| dt — i. Pe LAC CnO ae 
< [iv eo Aen 
i i Ipe(le — tI)I[{K-(e,t) — Ke(y,t)}ur())nl at 


forn > 1, r=1,2, we get 


b 
/ bp (lor — tl) — pe(ly — t))| [Kr (a, t)ttr(t)]nl dt << Mx (p,,h), 


b 
/ Ipr(la — e))|[{ Ao (a, t) — Ke (y, t)}ur(t)|n| dt < max |Ky(y,t) — K(x, ¢)] ||prllefa.o)- 


a<t<b 


From the uniform continuity of K,, for r= 1,2. With the same process, we can prove that 
|Torn Ur (x) > Tor ntr(y)| < é. 


Then, we get that the sets 


{Tirnt Ur € X, || er | |x < 1, nm = i} ? 
1<Ur<2 


are equicontinuous. 
Lemma 3.3.1. For n> 1, the set {anon > i} is collectively compact. 
Proof. We show that the closure of 


S= {An.0, n>1, Ux, < if, 


is compact. It is easy to show that 


S= TH Sten n>, ||| |x < i} ’ 
lr 


l=1 \r=1 


and by using theorem 3.3.1 and the Arzela-Ascoli theorem, we deduce that 


2 
{PTinntin n>1, |lurllx < i} 


p=1 lr 


is compact. Also, we have 


2 2 


2 
T]{SoTinner nr > Li, | ee | |x < a7 =[1{Soti-a nr = 1, \lur|lx < < i} 


l=1 Sr=1 l=1 \r= Ir 


Finally, the set S is a collectively compact. 


Theorem 3.3.2. Let Ar be a block operator matriz and Ar, be a block approximation operator. Then, 


VU € Xp, Jim |(Ar — Ar, UI |x, = 0. 
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Proof. For n large enough, we have 
b 
|Ty,U,(x) _ TirnUr(2)| < / \pr(x _ t)| |[A?(2, thur(t))n a K,(2, t)u,(t)| dt, 
such that for n >, r= 1,2 and 7 =0,...n 


[Kr (a, t)Ur(t)]n — Kr (a, t)ur(t) (ay —t) ie @j-1)Ur(#j-1) — Kr(z, out) 


? 


+ £(t— 25-1) ie 23) Ue( 3) — Ke(e,t)un(0| 


lA 


2 [Mero + |[Kco,0(Kr, lll 


On the other side, we have 


|Portir(2) — Tarnttr(a)| < i. Ip(« — £)| |[K (2, t)ur(t)]n — Kr(,t)ur(E)| dt 


b 
OK, OK, 
+ / Ipe(w— 1) [5 @, url) — F*(@,tur(6)| at 
such that 
OK, OK, 
[Fe (x, t)ur(t)Jn — is (z,t)ur(t)] < 2| Mao + ||Ko0,0(02Kr, h)| 00 el | ello 


+ 2| Maoh + ||Kico,0(Kr, h)|Lo0 | Ilp+-ILz3{0,0—a)> 


IA 


[AMtors ole h) + 2||Ko0,1 (Kr; h)|]oo jes I|Pr|lw2.2,0,e—a)- 


Finally, 


lA 


|(Ae= Ar, )U||x5 2 an (: he 2 lhwsoe-s) + M2 (: oF palhesoe-ci) Koo,1(u, h) 


+ 2 peak, Al (: + 2 lhwsoe~a) 


I|Fioo,1(K2, h)|]o0 (: zi 2palhwsoe-a) |U|lx,: 


When n — oo, we have h + 0 and we get the result. 


From lemma 3.3.1 and theorem 3.3.2, we obtain that our block approximation matrix verifies all of 
properties (i)-(iii), then by lemma 3.3.1, we get 


\||(Ar, — Ar)Ar||| 70, ||\(An, — Ar)Az,||| > 0, when n > +00. 


Theorem 3.3.3. Assume that 


Jim ||(Ar — Ar, Ullx, = 0, 


and 


\|\(Ar, — Ar)Arl||| 30, |[|\(Ar, - Ar)Az,||| 70, whenn- +00. 
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Then, for sufficiently large n, the approximate inverse (AIz — Ar,)~! exists and is uniformly bounded. 
We also have 


: 1+ [Nlfe— Aryl Aral 
Alg — A NiWh< = <0Oo. 
HOA Ar) IS KAR — Ar) MM (Ar — Ar, ) Anil 


Proof. See theorem 11.4.4 in [5]. 


Theorem 3.3.4. Let U = (wu, uz) be the solution of the system (3.4) and Uy, = (tin, U2,n) ts the solution 
of (3.5). Then, 


IU — Unllx, = 0. 


Proof. For n large enough, we have 


MU -U,) = ArU — Ax Un, 


= ArU Ar, U Ar, U Ar,, Un, 


= (Ar _ Ar, )U + Ar, (U = Un). 
Then, 
U —U, = (Aly _ Ar,)~! (ax _ A, 0] F 


which gives 


IF — Vallx, SMNOAL — Ar.) "Ill (Ar — Ax, UI Ix,5 


and when n — oo, we get the result. 


3.4 Collocation cubic b-spline method 


The purpose of creating a new method for solving this equation is not only to improve the approximate 
solution but also to try to build a fast method in solving it. In the latter method, we obtain an algebraic 
system of size 2n + 2 equations and when n is large enough then the system size will be large and thus 
the solution spends more time in computing. The tool we use to solve this problem is to apply one of 
the projection methods but with a choice of differentiable and continuous basis. So we take the cubic 
b-splines which are functions of class C?. The advantage of this method is that it eliminates the number 
of equations by half. Instead of solving a system of 2n + 2 equations, we solve a system of n+ 1 algebraic 
equations. 

So we apply the collocation processes and we need to give the definition of cubic b-spline functions (see 
[7, 18, 22, 32, 44)). 


Definition 3.4.1. Let A, = {a = 2 < 4 < ...% = b} be a uniform partition of interval [a,b], with 
Ra 
vu, =atith andh= ee B3(A,,) be the space of cubic b-spline functions, given as follow: 
n 


B3(A,) = {s € C? [a, b S |fesaeyi]€ P?,i=0,1,...n— i}, 


where, S |tx; 0:41] 18 the restriction of the spline function S : [0,1] —> R in each sub-interval [x;, 7:41] 


and P3 is the space of cubic polynomials. For i = —1,0,...n,n+1, we give the uniform cubic b-spline 
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from (4.5) 
(Go y9)*, OS [tea 874l, 
(x — x42)? — 4(@ — x4_1)3, xe [a;-1, v;], 
BS(0) = 75) (tne 2) ACen — 2),  € [i241] 
(vise — 2), x E [Li41, Vi42], 
0, otherwise. 


It is clear that {B_1, Bo, Bi,...,Bn—1, Bn, Bnsi} forms a basis of B3(A,). Then B3(A,) ts a finite 
dimensional linear sub-space of C?[a,b] with dimension n +3. 


We define (P?)n,cn~ the sequence of linear projection operator 


P?:C{a,b] —> Bs(An) 
n+1 
v r+ P3y(x) = x a; B? (2). (3.6) 


i=-1 
It satisfies the following interpolation condition for all v € C'[a, bj 


Peo) Sule). FH O01... 


(Prv)'(b) = v'(d). 


We construct an approximation solution u, = P?un, verifying the following equation 


Nn = P? Tun + P3f. (3.7) 


3.4.1 Convergence Analysis 


To study the convergence of the method, we recall what is spline interpolation. We will go through 
steps, including the passage from writing as a linear combination to writing using spline interpolation. The 
proof that converge to the identity operator, and that our approximate operator converges in norm.Finally 
the use of the latter to prove that the exact solution converges in the space C[a, b]. 


For the first time, it will be demonstrated that P?u convergence to u in the sense of the norm C'{a, b] 
especially since a lot of researches like [8, 49, 63],[33, 52] and [74] have shown this convergence, provided 
that is u is defined in C*[a,b], u € C°[-1,1] [80, 4] or it is necessary that u is periodic function [6]. 
Moreover, to prove this convergence, we recall the spline interpolation of u on the grid A,,n > 1 which 
denoted by S;,(x) (see [32, 33, 49, 52]) and defined as: In each interval [x;-1, xi] the spline S,,(7) has this 
formula which based on the moments M; = $”(a;) 


ge Mi (85) 
+ (ues mE) (4). (3.8) 


However, any cubic spline S,,(x) constructed on segment [a,b] is described on the linear combination of 
cubic b-spline [61, ?, 80]. Then, S,,(z) has the following form 


ey =o 


Va € [a,b], Sp(x) = Pe u(z). 
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This implies that: 
(ZL — Pryullerja,o = [lu — Srllerja,o- 


We return to the spline form (3.8) and for the left and right limit of the point x; 


_ h h u(a;i) — u(a;-1) 
' (x; = —M;,- M; 
Si (2; ) 6 1 + 3 + h ’ 
h h u(x; 1) = u(x;) 
S(at) = —--M,;-——M; ee 
rn (2, ) 3 6 414+ h 


The continuity of S/ (a) at x; yields for i=1,2,...n—1 


M;-1+4M; + Mist = (se) _ u(ai) u(x) = eo) ) 


h? h? 


(3.9) 


But, in many applications, it is more convenient to work with the slopes m; = $"(a;) rather than the 


moments M;. Here, we present another formula of S,, in each segment [a;_1, xj] 


Sn(z) = mia (x; — x) a — £i-1) ii (x — mea) (a; — 2) 
fe SE 
+ u(zi) (x — 24-1) etn —a)+h] 
Se) = ames (Ce = + a — 32) = (a — Nes ti-1 — 32) 


- 2 4 u(ax;) —_ u(xj-1) 
Sa, oS pi-1 + pMi 6 7 ; 
—4 2 U(ej41 — Ula; 

S"(e7) = ya Mitt 6 ( See (2%) 


For results, we require 


UX; — UX; 
mina + Ame + mig, = 3A elt at ee) 


with the end conditions 


Then, system (3.9) has the form CM = D and system (3.11) has the form Cm = E, where 
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(3.11) 
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m = (mo,m),... 


d; are given by 


2 1 0 
1 4 1 0 0 
0 1 4 1 0 0 
C= 
OO Ase SSa> wea OL 4 1 
2 se 25h “ae hee WE 22 


h ? 
ules) —ulas1) | ,U(®i41) — u(2s) 
& 23 “4 3 5 . $S 0 Het, 
a 3 U(tn-1) —u(arn) 


h ? 


and the vector E = (eo, €1,---€n)’ € R"*! have the following elements 


— ; (ae - u(Xo) wa) 


e; =6( Men ums). Aa) a), at a ee 
ce = 8 oy) Med =tne) 


Theorem 3.4.1. Let Pu defined by (3.6) and u the exact solution of (3.1). Then, 


I|u _ Peullora] <e Ki(u, h), 


where c is a positive constant. 


Proof. On each interval [x;_1,x;] and we use the spline form (3.8), we get 
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P? u(x) — u(x) 


I 
D 
3 

a 
8 
& 
| 
iS 
= 

8 

& 


6h 6h 6 
a, Faifgp |) | Po er) _ aes 
v 6 h ? 
— Mii 3. Mi-i,2 Mi 3 Mii>o 
aa (ai — x) Bh h’ (a, — a) + Bh (a — a; —1) 6h? (a — @-1) 


’ u(xi) + u(ai-1) aj+aj4-1—-—22 
| ( 5 fte)) (wea = u(r) | a 


We have 
(x; -—a)(a—ajz_1) < hi’, 
(2a;-—a;-1-2) < 2h, 
(a@—2x;-1+2;) < 2h, 
|a;taj1—-2¢| < A 
Then, we get 
|Sa(2)— ule) <P] Mal + ae f+ | SEHD Gey 


af | u(xi) fli) f 


this implies 


2h? 3 
Sn —Ulloo < =—||M|[pnti + = Ko(u, A). 
3 2 
On the other side, we have M = C~!D, then ||M||pn+1 <||| C7? ||| ||D]|pnti, where 


|Dllawrs = amax | di | 


3 
< he Ko(u, h). 
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Finally, 


“ 3 
Sn — us < (217 I +3) xolus 


Now, from the second interpolation formula (3.10) in each segment [2;-1,2;], we have 


n u(ai) = u(xi-1) _ ; (xi = x) (224-1 + 43- 32) ; (a = Ui-1) (24; + 4i- 3x 
Bal) h ae he ve re 
F u(xi) = u(xi-1) ‘ ; ; u(xi) = u(xi-1) 
t hh 6(x Li-1) (Xi x) h y 
with, 
6 @L— Vi-1 1 
6(a—a-1)(ui-—2) = h2 ( 5 ) 9? 
ee) t1t+ x Ld Cit ei 
(a; — ©)(2a;-1 + % — 3x) hp (: 5 ) a ( 5 ) 
So that 
S! (a) u(x;) ain) = E ( “J *) | mina — ul (ai_1) 
+ (m;—1u'(a;) -2 (“4 ain) )| 
: (« Li “ct) o tte 
— w(aj_1tu' (a) + w(a)| : 
Then, 


| (x) — Mm) IY | mia — ul int) | + | me — (as) | 42 ou! A) 
+ |m,—u'(a;) | + | mii — u'(ai_-1) | +Ko(u’, h). (3.12) 


We denote U = (u(xo), u(x1),...u(ap,)) and U, = (u’(xo), u'(#1),...u’(a@,)). In other way, we have 


1 E 
rm — Wages SI] = EE |ts | SEO Uo (3.13) 


Also, we have 
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where, I;,41 is the identity matrix of size (n +1) x (n+1). We have 


6 us 3 6 
En-1 En 
6° 8 
with the end and the first conditions, we get 
| oe <3 Ko(u’,h), 


| —Ej-1 ej Ej41 | 
6 3 6 


<6 Ko(w',h), i=1,2,... 


TEn— En 
Igo t gl <8 xo(u',h), 
which gives || E|pn+1 <6 Ko(u’,h). Then, from (3.13) 


E = 
rm = Slants <6 ||] C77 [I] wo(u’ 


Substituting (3.14) in (3.12) 


| She) — SEO MI) fe (12 I] o7 I] 41) solu 


and therefore, 


We get 


/ / / u(a;) = u(a—1) i , 
[Sa —w'lc << max, | Sale) ee) T+ mag, | Me u'(a) | 
< (u |] Oot III +2) Ko(u’, h). 
Then, we obtain from (3.12) and (3.15): 
Iu _ P2ullora] = Iu _ Snllor[a,b] <e K1(u, h), 


3 
where c = max (11 ||| C2 ||] +1, 2 ||| C7? ||| +3). 


(3.14) 


(3.15) 


Theorem 3.4.2. Let P? be a projection operator given by (3.6) and T is a compact operator defined by 


(3.3). Then, 


, _ ps = 
lim ||(I— P8)T|| = 0. 
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8.4. Collocation cubic b-spline method 


Proof. Since T is compact the set 
M = {Tu € C’|a, b, Il24I| 4 [a,] < i} 


is relatively compact in the Banach space C*{a,b] and by Banach-Steinhaus theorem [5], P? converges 
uniformly to the identity operator I in C'{a, | 


lim ||(I-P;)T||= lim sup || — Pz) Tullorfa,s) = sup || — Pr)vllorta,s}- 
n—0o n> veEM 


a ell ot fa.) 


Under above theorems, we get P? is pointwise convergent to the identity operator and P?T convergence 
to T. Then, (AI — P3T)~! exists and ||(AI — P3T)~+|| < 00 (see [81]). 
Theorem 3.4.3. Let un, be the solution of (3.12) and wu be the exact solution of (3.1). Then, 
lim ||u — un||crfa,o] = 0 
N—- oo 


Proof. For n enough large, u— un, = (AI — P3T)~} jir- P3T)u+(f- P31)| and we obtain the following 


estimation 


Iu — unllorfa,) < AL — PT)" | je Pee le= PS) lees) 


When n — ov, we get the result. 


3.4.2 System Approximation 


From the equation (3.12), we have for all x € [a, }] 


un(2) = P?Tun(x) + P3 f(x). 


We put x = x; for j = 0,1,...n, then by the conditions interpolation of the sequence (P?)nen (3.4), we get 


Aun(zj) = Tun(aj) + f(x;), 


which is equivalent to the following algebraic system for 7 = 0,1,...n 


n+1 b b 
ya [2) — f pulls ~ t)Ka(es,t) R(t) dt — f pal| ny — t) Kaley, BS (bat 
i=—-1 a a 
=  f(zj), (3.17) 
where {a,;}"*!, are unknowns coefficients to be determined. 


Study of linear integral differential equations of the second kind: Analytical and numerical approach 66 


3.5. Numerical tests 


But we have a problem with the system (3.17) which has n +1 equations with n +3 unknowns q;. For 
this reason, we introduce B?;(x) be a modified Cubic b-splines basis [31] 


B3,(z)= B3(x) +2B%,(z), for 7 = 0, 

B3,(z)= BP (a) — B3, (a), for j = 1, 

Bi (¢)= “B(a), for j = 2,2... — 2, 
BS, 1(¢)= B3_,(x) — BR, (2), ior 7 == 1, 

B3,(2)= B3(x)+2B2,,(2), for j =n. 


In the next table, we present the value of B3;(2;) and its derivatives 


Table 3.1: Value of modified b-cubic spline and their derivatives 


x; j=0 jfzi-2 jsi-1l j=i jsitl j=i4+2 jeHen 


B3,(x;) 6 0 1 4 1 0 6 
—3 3 6 


Finally, the solution u,, has a new representation 


Va € [a,b], un(x) = ai B, (x), 


and the system (3.17) has a new form given by 


b b oy 
y a;] ABi°(e) - [P rvlles =) (2j, BP (0 at — [pj —t |) Kaley 1) BE at 
= f(zxj), 


where {@;}"_) are a new coefficients to be deteminated. 


3.5 Numerical tests 


To illustrate the efficiency of our numerical processes and to compare between the both methods, we 
give two numerical examples. The first is defined as 


3.5.1 Test O1 


Consider the following integro-differential equation 


u(2) = f Vea art f |x — t|Fu'(t) dt + f(x), (3.18) 
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where, 
Qr+12)(1—a)? 8a%e7 9rd | 3 8a? + 124 +15 
= 7x? ( al 2 7 
5 al 14 7 ae ge 105 
We have 
P,(je-t]) = Vle2—-dl, 
P,(|x-t!) = |e-t]3, 
and 
ev 
Ki(z,t) = 
1(2, ) t+? 
Kx(\c—t|) = 1. 


Moreover, M, = e and M2 = 1. In addition, 7 > e*||Pi||w+.1[0,o—a) + ||P2||w2110,0-a), then the equation 
(4.9) has a unique solution. 


The next table present the error discrete between the exact and approximate solution which given by 


Try = Bax |un(wi) — u(ax;)|. 


Table 3.2: The error between the exact and approximation solution of equation (4.9) 


n Product integration time(seconds) Collocation cubic b-spline time(seconds) 
10 2.4205e-04 0.985059 4.7601e-05 0.620084 
50 9.6603e-06 14.386013 8.5975e-07 7.456484 
100 2.4147e-06 51.016760 1 .6386e-07 27.258581 


Study of linear integral differential equations of the second kind: Analytical and numerical approach 


68 


3.5. Numerical tests 


x104 


2.5 1 


—— B-spline Error A 


Product Error 


Figure 3.1: Errors of B-spline and Product integration method (4.9) n=10. 


3.5.2 Test 02 


Consider the following integro-differential equation 


10 10 
u(a)= fo fe-tisut) at+ fo Vie tulle) at + F(a), ie 
0 0 
where, 
‘i 3 2 3 9 ee 
f(a) = 2 — = (8x - 10)(10 - 2)$ — (10-2) + Sah + al 
We have 
Pi(le-t)) = |e—tl8, 
P,(le-t]) = Vie-t), 
and 
Ky (z,t) = 1, 
K2(|x — t}) = 1. 


Moreover, M, = 1 and Mz = 1. In addition, 11 > ||Pi||w2.1[0,o—a) + ||Pallw210,o—a), then the equation 
(4.9) has a unique solution. 


Table 3.3: The error between the exact and approximation solution of equation (4.4) 


n Product integration time(seconds) Collocation cubic b-spline — time(seconds) 


50 1.9077e-07 14.038054 5.9756e-11 7.198711 
100 1.8224¢e-07 57.896836 2.2340e-11 33.793924 
200 4.5985e-08 213.576772 8.7540e-12 116.348824 
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T 
—— B-spline Error Y 


1.87 — & — Product Error 


/ 


“Oa ppqqcoEoS? Be, is A 
aie "SEBOBO GG. A 7 


Figure 3.2: Error of B-spline and Product integration method (4.4) n=50. 
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Chapter 4 7 


Analytical and Numerical Study 


of Linear Fredholm 
Integro-differential Equation 


In Sobolev Spaces 
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In this chapter, we study the existence and uniqueness of the solution of a linear integro-differential 
Fredholm equation (2.1) in the most complex and weak Banach spaces, which are Sobolev spaces. We 
suppose a theory that proves the existence of the solution under certain conditions, which allow us later 
to be in harmony with our two analytical and numerical studies. This theory is well detailed in section 
4.1. In this chapter, we are going to construct two solutions, each based on two Projection methods: 
Galerkin and Kantorovich which are well explained in section 4.2. We present some theorems through 
which we will prove the convergence of both solutions to the exact solution, but only in the sense of the 
norm of two spaces H'{a,b| and W1![a,b]. This what we present in sections 4.3 and 4.4. Finally, we 
give a numerical examples 4.5 to see the error computation of the Galerkin and Knatrovich method in 
the two Sobolev spaces. 


Let p € [1,-+co|. We recall the linear Fredholm integro-differential equation 


b b 
Va € [a,b], Au(x) = / Kyi(a,t)u(t) dt +f Ko(a,t)u'(t) dt + f(x), (4.1) 
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in the Sobolev space W'? is equipped with the following norm 


lol |v. pao) = [elle fae) + [Yl ce fa,8)- (4.2) 


4.1 Analytical study in the Sobolev space W!?{a, }] 


In this section, we are interested on the solution existence and uniqueness of (3.1) in the Sobolev space 
W!?/a,b]. We propose two different hypotheses on the kernels Kj, i = 1,2 depending on p € [1, +00 


(£1) For p €]1, co[, we assume that the kernels K;, 1 = 1,2 checked the next assumptions 


(L1) | K;(a,t) € L?((a, b]?,R) and (01 € L([a, b]?,R), 


(£2) For p= 1 suppose that 


b 
| ames [1K |i; (a, t)| dt < +00 and max, [ 
a<a<b <b 


OK; 
Da 0) dt < +00. 


Now, for all p € [1,-+oo[ and u € W! [a,b], we define the following linear operator 


Va € [a, b], a= ff Ky (a, t)u @ ars [ Kele,OW( dt. (4.3) 


Proposition 4.1.1. Let p € [1,+o00[, f € W! [a,b] and K;, i = 1,2 verify one of hypotheses (Li) or 
(Lo). Then, Tu € W'?[a,b] for allu€ W'”. In addition, the weak derivatives (Tu)' is given by 


K 
a= [Ue *(x,t)u oar fe co *(x,t)u'(t) dt, p.p for all x € {a,b}. 


Proof. To prove Tu € W? [a,b], it is necessary show that Tu € L?{a,b] and (Tu)’ € L?[a, bj, for 
p & [1, +00. 
First, we verify that Tu € L?{a,b] for all p € [1, +00] 


I. For p €]1, +00|, we use the Hélder’s Inequality [1]to deduce 


[ise.ow piace (fo IK(e.0le a)’ ([ lu(t ye at) (4.4) 
[satu nate (fi ale,0)¢ a)" ([ lu (t yr at) (4.5) 


By combining equalities (4.4) and (4.5), 


| [ Kate.nue aur [ Kaledu'o at |< (ice, a) ¢; ju(t yr at) 
+(/ "| Kale, 0) at) (/ Wor it) "46 
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Then, 


Pp 


< Cxcor it) ([ |u(t) |? a)’ 
+([ inateo w)'(f ju! (e)|P w)'T 


he Ky (a, t)u(t) dt + [ Ko(x, t)u'(t) dt 


Applying the previous results, we get 


b b a b a7P 
[Tell ofa,0) < | I(/ Jka (e.0) Pat) F (/ IKs(e, 0) | dex ||u|| Ww. [a,0)- 
Using the fact (z+ y)? < 2?-1(z? + y?), we have 
Tulle) < 1) ([ |Ki(z, t)|? a)’ dx 
+ [ ([ Kaeo at)" te || |we2.2fa,b] < +00. 


x 


II. For p= 1, we have 


b b 
[IPullertae < (nas | [Ki(x,t)| de-+ max, / |Ko(c, 1) it) [Pull wa.sfa,s) < +00. 


a<a<b 


This proves that Tu € L?[a, b], for p € [1, +00. 
Secondly, we show that there exists a function g € L?[a, bl, p € [1, +00] verifying 


b b 
; Puede) da = / g(x)4(a) de, Vo € C™[a, 8]. 


Let ¢ € C™[a, b]. Then, 
[tune os [[ [ seoue ar [Kaloo | boi 


= [[f seoe 1) de] u(t as fo | Kate. 90% 2) de] (t) dt, 
— [| f Beno ae] nae le ™ (0, 1)d (2) de} u ‘(t) dt, 
= -[[f (0, tut nay [Xe as * (x, t)u (0 at| (0) de, 


7 ii (Tuy! (a) h(a) de. 


I 


Hence, V¢ € Coa, b], we get g = (Tu)’. Finally, let verify that (Tu)’ € L[a, b). 


q PR 
qd 
it) ae ||| wasp [a,b]: 


I. For p €]1,+00|, with the same last processes. We can prove that 


1era'ttnoa sf f (f° a) ae f'(f 


Ok, 


OK2 
Se (a;t) 


Ox (24) 
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II. For p= 1, we have 


b 
/ 
[u)'Ihesto S (soa, i: 


So, we have Tu € W![a, b], for all p € [1, +00], this completes the proof. 


Oki 
Ox 


dt + max 
a<a<b 


volar) |e] | w2-1[a,b)- 


(x, 


Under proposition 4.1.1, the derivative u’ is given implicitly as 


Vae€[a,b], Au'( a= [% ies (x, t)u oar fi 2 (0, tu (0 dt + f'(z). 


We introduce the following theorem to show that T € BL(W'?[a, J). 


Theorem 4.1.1. 


1. If the kernels K;, for i= 1,2 verify the assumption (Li) and p €j1,+o0[, then T € BL(W1[a, b]) 


and 


uy 


7 \qG yp 
at) ac| ; 
then T € BL(W'[a, b]) and 


eee (2.0 at), 
x 


OK; 
Ox (a, t) 


nisl (f wdeor we)’ w+ f(f 


2. If the kernels K;, fori =1,2 verify the assumption (Lo) , 


2 b b 
< ; 
\rI|< > (2, [ isie.olae+ amas, | 
i= 


Proof. We recall that the norm of linear operator T is 


IZI}= sup | [Tull w+»[a,0)- 


Hel wt.p [ao] = 


(4.8) 


1. We start by prove that the operator T is bounded in W [a,b], for p €]1,-+oo|. By summation of 


(4.7) and (4.7), 


2 b b 
[Pulleoray FIT) Bray < a 3 / ( Fi 1Ki(o.0) at) Fe 
i=0 . sad 
b b Bg 
“LU 


x 


Ky 
ao (w, 8) 


2 q 
it) ae ||ul|w1.2[a,6)- 


On the other side, we have 


p 
(irreliren 7 (TeV lion) <ort (iret, y+ Me) eo fa, 4) 


We get, 


oe b b c 
>| / (/ Kile, at) dx 
+ Le 


Ox 
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||Tul|z>ja,e) + (Tu) '||22f0,0) 


(x,t) 


q 2 A 
ar) ac| ||| |W. [a,b]: 
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Finally, 


2 


ge b b a b b 
\Tullwiopey $222 >| / (/ IKi(c, 9) dt) ax+ f (/ 


i=0 


OK; z 
Het) iv)’ ar)’ fellas etaa 


This last equality proves that Tu € BL(W1?[a, b]) for p €]1,+00[. Then, 


2 


Z| <9) f( fixiane we)’ w+ f(f 


i=0 


OK; 
Ox (x, t) 


co \? 4h 
ar) ae|”. 


2. Now, we prove that T € BL(W1[a, d}). 
From (4.7) and (4.10), we find that ||T'u||yw11,,9) is increased by 


2 


||Pul|w.1[a,0) < (ms max x, [Us | A; (ax,t)| dt + max, 


i=0 


co) | ae) ||| | w72-2 [a,b]; 


so that Tu is a linear and bounded operator in W![a, b]. Using the norm of linear operator, we 
get 


2 b b OK, 
< 4 
TI < > (295, [ise piace max, [| (2.0 at). 


In the next corollary, we give a sufficient condition to ensure the solution existence and uniqueness of 
equation (4.1). 


Corollary 4.1.1. Equation (4.1) has a unique solution in Wt? [a,b if 


1. For p €j1,+col, the kernels K; for i= 1,2 check the assumption (£1) and 


ELL (Losers) as ff 


1=0 


OK; 
5 (x,t) 


ou 1 \G de 
|A| > 2°? it) as| (4.9) 


2. For p=1, the kernels K; fori =1,2 verify the assumption (L2) and 


2 b 
al> O(a, [ |; (a, t)| dt + mua, [ 


(x, 0] at) (4.10) 


Proof. 


1. To prove (4.9), we will show that 


LUE Limeara)'on LC 


i=0 


OK; 
5 (x,t) 


2p— 
[T|| = 2°» 


a \t 9h 
it) ae”. 
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By using the sup definition, we have 


Ve>0,du1 €W'(Q); |lur||wrrfa,s) = 1, ||[Luallwrfae] > |||] -e, 


since, 


1 


q q 7 
ir) as| ||u1||w-»[a,b] 


(x,t) 


$f Lmeora) f(s 


x 
> ||All—<«, 


For any arbitrary ¢, we get 


Ef (fuscnra) aes f( 


OK; 
t 
Ox (z, ) 


q PR 1 
at) as| > ||T|I, 


and 


2 


STL (fixer a)! aes ff 


i=0 


ied : 5G Vs 
\|T\| =2°7 5 eu) it) as| < |Al. 


By Neumann’s theorem 1.7 (AJ — T)~+ exists and 


1 


[OLS 2" I <5 ea 
|Al — II7| 


where J is the identity operator defined in W1 [a,b] in itself. Then our equation has a unique 
solution. 


2. To show the equality (4.10), we will prove that 


2 


b b 
I= > (2s, [ise ol are imax, f 
i= 


Also, by sup characterisation we have 


OK; 
= fu) at). 


Ve>0,du €W(Q); |lurllwrfas) = 1, ||[Puallw.1fa,) > |All -«, 


such that 


2 b b OK; 
(may, f 1Ki(at)| db+ mass | plat} at) Telnaes 


> ||Tur|lw2rfa,s > ||P I] - . (4.11) 


Since ¢ is arbitrary and from the equality (4.11), we get 


b 
WPM = (25, [ |; (a, t)| dt+ max 
ax<a<b 


OK; 
Da (=. it) < |Al. 


Using Neumann’s theorem 1.7, we get (AJ — T)~1 exists and bounded, such that this time I is the 
identity operator of W!{a,b]. Then, the equation (4.1) has a unique solution in W!1{a, b. 
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In the rest of our thesis, we present different constructions of numerical solutions and we show that 
these approaches converge to the exact solution. So, in the convergence analysis we need the compactness 
of the operator A, which is proved in the coming theorem. 


Theorem 4.1.2. Let p € [1,+co|. Then, T is a compact operator. 


Proof. Let {un }nen be a sequence defined in W!?[a, b], since ||un||w2.r[a,b) < 00. {T'Un}nen is bounded in 
W'[a,b]. Therefore, it has a sub sequences {Tun, }gen converges in W1? [a,b]. Consequently {Tun }ne 
is precompact. Thus, T is compact. 


4.2 Projection Approximation 


In this section, we present the basic idea of projection methods 1.4.1 in the Sobolev space W1?{a, 0]. 
We put X = L? [a,b] and we denote Xz by X x X. Let {Xn}n>1 be a sequence of finite dimensional sub- 
space of L? a, b] and {Xon}n>1 be a sequence of finite dimensional subspace of X. Here, the construction 
of an approximation solution in the Sobolev space is based on the following projection operator {P,}n>1 
which is given by 


P., : Xo — Xan 
Vor> PrV a (Pi avis P2 V2), (4.12) 


where {Pyn}n>1, for r= 1,2 are sequences have the following form 


YWeEX, Prynv= So, e; ex, r= 1,2, 
i=0 


such that, (.,.) is the duality product in Banach space, {e;}7_) is an ordered basis of X,, and {e; }"_, are 
linear continuous functional, since 


To illustrate our numerical process in W1 [a,b] by using the Galerkin 1.4.1 and Kantorovich 1.4.1 
approximation methods, we need to present our operator T in the equivalent form. For this reason, we 
set u’-)) = u, for r = 1,2 and we define the following block operator matrix Ay which is given as 


Ar 7 re —- Xo 
U =(u1,u2) > ArpU = (Trim + Theva, To1ui + Th2u2), 


with the operators {Ti, }i<ir<2 are defined for all x € [a,b] and r = 1,2 by 


b 
VeeR HaG). = i K,(a,t)u(t) dt, 
> AK, 
Vue X, Teyv(x) = 5 (a, t)u(t) dt. 
a vo 


We rewrite equations (4.1) and (4.8) in this simple version: 
(Ala _ Ar)U = F, 


where Ip is the identity operator of X> and F = (fi, fo). 
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Define now, two matrix approximations AF = = P, ArP, and AK = P,Ar. The first is called block 
Galerkin approximation matrix and the said one is block Wontoevicl matrix. Using the above ap- 
proximation operators Galerkin and Kantorovich, we get two approximation equations. 


1. The first one is the Galerkin equation, which has the following form 
(\Iz — AG Un = F. (4.13) 
From the equation (4.13), we have 


I 


PB Ap Pi ntins£2,,tan) + (Fis to), 
= Pi(TuPintajyn + Ti2P2ntvan,T21Pintijyn + T22P2ntan) + (fi, fe) 
= (PinTi Pinan + PonTi2P2ntv2,n,PinToiPintijyn + P2nT22P2nU2,n) + (fi, fa). 


Thus, the unique solution Up, = (ui1,n, U2,n) of Galekin equation (4.13) is given by 


[Ed Came 2) (Tas Bae Y Gage) Ones, ce tA. 
~ a 0 1=07=0 


U2n = 7 5 [Ed ume e5) (Taie;,e eres + 32 (tans e5) (Tr2e;,e ee +e) 


1=07=0 


such that, Xi (i) = (vin,e%) and X2(i) = (v2, e%) are unknowns that we will find by solving the 
following algebraic system 


Xy Ay Ai Xy Y1 
At |= — |4t_ |, (4.14) 
Xo» Ag, Ago X2 Yo 
with 
Anir(t,j) = (Treief),1<ir<2, O<tj<n, 
Yi@) = (fe), DSc<n, 
Yo(t) = (fo,e7), O<i<n 


2. The second one is Kantorovich approximation 
(\Iz — AKU, = F, (4.15) 


From the equation (4.15), we have 


Mn P, ApU + F, 
= Pi(Tiitijyn + Ti22n, To1Uin + Te2t2,n) + (fi, fa), 


= (PinTiayn + PenTi2ten, PinToiuiyn + PonTo2tan) + (fi, fe). 
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Then, the unique solution of Kantorovich system has the following form 


1 n n 
Un = X b (Tia tig €; ec S Tietle 8; ee + f| ; 
i=0 


i=0 


hry = 5 (Stam, nee eet Y (Catt, ny € ete + fa] 
The unknowns of the previous system are Xy(i) = (Tii4in + Ti2t2,n, e7) and X2(t) = (Taitijn + 
T22U2,n,e;) for i = 0,1,...n. To find them we multiply the system by the block matrix Ay and 
apply ej [48]. We find this n+-1-dimensional linear system with the two unknown vectors X; and X2 


Xy Anji An,12 Xy Yi 
‘ a = (4.16) 
Xo Ani An,22 Xo Yo 
where 
An tr(t, J) = (Tines, €}), 1< Lf < 2, 0 < 14] < n, 
Yili) = (Tuf+Taf',e7), O<i<n, 


Y2(#) (Toif + Toof’,es), O<i<n. 


I 


4.3. Study of Galerkin and Kantorovich approximation in W'[a, }] 


In this section, we set X = L'[a,b] and X2 = X x X. Let n> 1, we recall A, the uniform partition 
of [a, 6] (1.15).Let Xp.2 = Xn X Xn, where X;, be a subspace of Dia. b] has a basis {e;}"_, such that 


1 if GE [vi-1, Xi], 
e;(x) = 


We define the duality product as 


Lemma 4.3.1. Let P,, : X. > Xon be a sequence of projection operator. Then, 
1. Py € BL(X2, Xan), 
2. For all U € Xo, we have 


(Za = Pall, $2 M1a(h,w) > 0, 


which means that (Pn)n>1 is pointwise convergent to Ig. 


Proof. It is clear that {P,}n>1 is linear. We need to prove that it is also bounded. Since 


1 n b Li 

5 ff butailes(a)| dy ae, 
t=O % “M1 

ieee Ly Ly 

om jus (y)| dy da < (b—a)||ual|x, 
i=0 % Vi-1 M Vi-1 
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I|PnU||x, < (6—@)||U||x,- This implies that {Pp }n>1 is bounded and ||P,|| < (b— a). 


We have, for n > 1 


[|U — PU ||x, = |] — Pin)urllx + [|Z — Pen)uallx- (4.17) 
But, 
b 
W(—Panduallx = fh lua(0) ~ Pata()| dr 
= / |ur(x) — Py nus(2x)| dx 
i=1° 7-1 
T= Li Li 
= wf / |us (x) — uy (t)| dt dex, 
i= Ti-1 VY Ui-1 
< Sy "Jus (2) — uy (t)| dx at 
U4\ Lv) — xv 
SS he ben) 1 U1 , 
9 n Li h 
< </ finale +8) ~wa(a)| ae at, 
i= Li-1 0 
9 b h 
< 5 ff e+9-ul@)|arae 
h a 0 
9 h b 
<5 ff iue+t)~m(o) dtae, 
h 0 a 
< 2 K1,0(A, u1). (4.18) 
In the same way, 
\|(I _ P2n)ual| x < 2 K1,0(h, ua). (4.19) 


Substituting (4.18) and (4.19) in (4.17), we obtain 


\|U — P,U |x, <4 Ki1(h, U), 


b— 
so that for h = pee when n — +00 we get h > 0 and 
n 


lim ||U — PrUllz, = 0. 


4.3.1 Kantorovich Approximaion 
In this section, we study the convergence of Kantorovich in the space W!1{a, b]. 


Theorem 4.3.1. Let Ar be a block compact operator matriz and Ak be a block Kantorovich approxi- 
mation. Then, 


(Ar — AX Uz, < 2max (Ish, Khe F104, Kalle) 


and we have lim \||(Ar — AF )IIl =0. 
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Proof. By using the norm of block operator matrix and the equality (4.17) in lemma.4.3.1, we get for 
IIx, =1 


2 
(Ar - AKU, << max DOIN Tip  7Tie)urllx 
a 


Fes Sal o(h, TipUy), (4.20) 
where, for r= 1,2 
ao(h,Tirtr) = sup if [Tiptip (2 + y) — Tirttr(y)| dy, 
xEl[0,h 
< up ff iKst 2+ y,t) — Ky(y,)| lun(t)] dt dy 
xE[0,h 
< sup 7 sup |K,(a + y,t) — Ky(y,t)| dy |luell x, 
x€[0,h] Ja té[a,d] 
b 
< sup f sup [Ky(e+yt) — Ko(yst)] dy [ull 
tE[a,b] Ja xE[0,h] 
< sup &1,0(h, K,)(t). 
tE [a,b] 
In same way, we get 
K1,0(h, To-tr) < sup K1,0(h, 0, K;) (t), (4.21) 


tE [a,b] 


where, 0, K,. is a partial derivative of K;, respect to «x. 
Substitute (4.21) and (4.21) in (4.20) to get 


(Ar — AF, UL x, 


IA 


2 max (ath K1)(t), ®1,1(h, Ka)(t)), 


IA 


2( ahs Bh i fieecte K}llx 


Let us prove that lim |||(Ar — A% )||| = 0. Let Ar be a compact bloc matrix operator. The set 
n—>0o is 


{ArV, ||V||x, < 1}, is relatively compact. Using the Banach-Steinhaus theorem, the pointwise conver- 
gence of P,, to Iz is uniform and we obtain 


: K ‘ _ 
fim (Ar — Az, Ill = [IlZ2 — Pa) All| = Aree alee — Pa)V|lx, = 0. 


Theorem 4.3.2. Let u be the solution of (4.1) and let un, be the Kantorovich approximation. Then, 


lim ||w = Un||w2 [a,b] = 0. (4.22) 
n—0 


Proof. We have that A converge to Ar in the operator norm. Then, (AIz2 — Af,)~! exists and is 
bounded ( see [2, 5] ). We can write, for n large 


U —U, = (Alp — AF) (Can = AK )(uu!)), 
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we set C2 = |||(Al2 — AZ, )~+||| to obtain 


I|u —unllwarfas} = MU -Ualle, < Call(Ar — AF, )(4, 4 )ILx, 


IA 


2 Coy max (anc K1)|loo; Fa, Kall ‘ 
Then, 


lim || tn _ ullw21[a,b] =0. 
N—-Co 


Kantorovich System 


Now, we introduce the elements of the system (4.16) 


1 Ly aa 
Ai,(i, J) a 7 K,(a,t) dt dz, 1<i,j<n,r=1,2, 
j-1 JLj-1 


~ | [ OK, 
Wee odo, OS 


Yi(i) = ; | [ [Kier as [ KeoF00 a| dx, 1<i<n, 


Ly 
Li-1 


Aay (i, j) 


(a,t) dtdz, 1<i,j<n, r=1,2, 


Li b b 
Yo(i) = a | f (0, F(0) a+ | (0,10 a| dz,1<i<n. 


a 


4.3.2 Galerkin Approximation 


In this part, we present some theorems to show the convergence of the Galerkin solution. 


Theorem 4.3.3. Let Ar be a block compact operator matrix and Ag be a block Galerkin approximation. 
Then, 


Jim |||(Ar — AG, )Arll| =0, and tim |||(Ar — Ar) AG,|I| = 0. (4.23) 
Proof. For n large enough, we have 
(AZ, — Ar) Az, = (At, — Ar) AF, 
and 


(AG — Ar)A = (AK — Ap) AK + Ar(AK — Ar) 


IAF, — Ar) Ag, III 
IIIA, — Ar) Arll| 
By using theorem 4.3.1, we get 


NM(AZ, — Ar IM Pall? IAL, 
(I|Pall +1) (IAzlll IIl4z, — Arill- 


IN IA 


IIAz, — Ar) Ag, |I] < 2(b — a)? ||| Ar||| max (Iina(h, Ki) F1a(8, Kalle) 


II(AZ, — Ar) Arl|] <2 (6- a+ 1)|||Ar||| max (anc Ky )Iloo» (8, Kalle) 
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and when n — oo, we get the result and we have that (AJz — AG)~! exists and is bounded (see [2]). 


Theorem 4.3.4. Let u be the solution of (4.1) and let un be the Galerkin approximation. Then, 


lim ||u — Un||w11[a,b] = 0. 
n—0 


Proof. For n large enough, we have 


U — Uy = (Alz — AG)! as (P,U —U) + (AR — Av, 


| — Unllz, < c1 iI ||PnU — Ullx, + II(AR, - Ar\UIlz) 


By theorem 4.3.1 and lemmad4.3.1, we get 


= Call, <2 ©, [MAR M ,a(,u) + mar (aC, Ba), 1 (0, Ka)(0)) 
Then, 


lim || tn — ullw21[a,b] = 0. 
n—- co 


Galerkin System 


The elements of system (4.14) are presented as 


Air(i,j) = if Been dtdz, 1<i,j<n, r=1,2, 
Ao, (i,j) = Fed. a (at) dtdxz, 1<ij<n, r=1,2, 
ni = f° fleas, 1<i<n, 
Y2(i) = if Pi@jds, 1<7<n. 


4.4 Study of Galerkin and Kantorovich method in H'{a, D) 


Let X = L?{a,b], and X,, be a subspace of L?[a, b] has the following orthonormal basis 


2x —b-—a 


Vane la, DI, e;(z) = V2i+1 L,( b— 


), 


such that L;(t) is the Legendre polynomial defined on [—1, 1] by the following recursive formula 


Lo(t) =1, 
Ii(t) =t, 
2i+1 a 
Li t = = Lilt : Lie t). 
n@) = Sh) - Soha) 
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Since for all v € X and j =0,...n 
b 
(i, Bs) =| u(z)es(x) da, 


the subsequence operator 7,, be an orthogonal projection in X. Then, 


Lemma 4.4.1. Let P, :X > Kee be a sequence of orthogonal projection, then 
Jim ||(U2 — Pa)UIIg, = 0, 
i.e P, is pointwise convergent to Iz. 


Proof. Let P be a polynomial of degree n. Then, we write p, as a linear combination of Legendre Poly- 
nomials for any coefficients a; 


Ve €[-1,1], P(x) =) lajL;(a). 
j=0 


Let us prove that 


(uy a Pint, 9055) = SCai (ur — Pint, Li;). 
j=0 j=0 
we have 
(u1—Pinti,Lj) = (ui, L5) — So (Pints, Li) L;, 


ll 
° 


= (ui, 5) — (ui, £;)(L;, D3) = 0. 
This leads to 


n 
(u1 = Pintn, > aye;) = 0. 
j=0 


By Pythagorian theorem, we get 


n n 
[tan — Prt, > ag L5)\IF27-14) = lea — PrteallZ2¢aj + [Pinta — Sag LyllFep ay: 
j=0 j=0 


Then, 


||uy = Pr nuallt2¢ 1] < | 


n 2 
(u1 = Pints, Sa; L;) 
j=0 L?{-1,1] 


From the density of C°{[—1,1] in L?[-1,1], we get for any function u; € L?[—1,1] there is a function 
v € C°[-1, 1], such that ||u1—v||z2;-1,1) < =. In other way, by Weierstrass theorem for any v € C°[—1, 1], 


there is exist a polynomial P such that |v — P||coj-1,1) < 5 So, 


[lv — Pllzzt-ayy < V2||v — Plicof-1,14) S 


oR) 
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Finally, 


A 


nm 
[ler — PantéallZ2¢13) S lla — Panta, Sag LgllZoeaayll < |lus — ollcet—ayay + lle — vllfepay 


j=0 


IA 
om 


In the same way, 


|| we = Po nte||t2¢1,1] < E, 


which gives 


lim [|Z — Px)UIlxz, = 0. 


4.4.1 Kantorovich Approximation 


Theorem 4.4.1. If {Pn}n>1 is pointwise convergent to Iz in Xo, then 


: K _ 2 
(lim |AK — Ar|| =0. 


Proof. Fist, we prove that Ar is compact. Let U, be a sequence of X2, such that |U||z, <0. We have 


ATU lz, < CITI lz, 


Then, Ar is bounded, so that it has a convergent sub sequence {A7rUn}n>1, which proves that Ar is 


compact. 


Since Ar is compact, then M = {ArU, ||U||x, < 1} are relatively compact sets in the Banach space Xo 


and by the Banach Steinhaus theorem’s ( [5, 2]) {Pn}nen converges uniformly to Iz in M, i.e 


aia ~All = Ag aw LALO Ar 
Xo 


lim sup ||(J2 — Pr)4Z\Ix,- 
EM 


N07 


Theorem 4.4.2. Let u the exact solution and uy, is a Kantorovich solution. Then, 


lim ||u _ Un|| #4 [a,b] = 0. 


Proof. For n sufficiently large, U — Un = (Alp — Af )~1(Af, — Ar)U. Then, 


|| — Unlleafa,ey = NU —Unllx, < Call(Ar, — Ar)UI Ixy. 


By the last theorem, when n — oo, we get ||u — Un||H1[0,5) > 0. 
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Kantorovich System 


From the definition of the basis (L;)"_9, the elements of the system (4.16) are given by 


b pb 
Ai,-(t,j) = / K, (2; fe(aje,(t) dide, 1 tj <n, r=1,2, 


1 f° f° aK, 
Aa(ind) = 5 ff FE lededaer(t) at de, 1 <i n, r= 1,2, 


b 
a 


Yay = [[[ se@onoen) arf Kle0F00 at} eo dee LG, 


bp pb b 
Y2(t) = / | (0, (0) dt + (0,40 a| ei(z) dz, l<i<n. 


4.4.2 Galerkin Approximation 


Theorem 4.4.3. Let Ar be a block compact operator matriz, Ag and Ag are the block Kantorovich 
and Galerikn approximation matrix respectively. If lim |||A4 — Al|| =0. Then, 
noo mm 


lim ||(A§, — Ar)Ar|||=0 and lim ||(A¥, — Ar) AG, || = 0. 


. Vv 
This means AG v-convergence ( AG —> Ar, ). 


Proof. For n large enough, we have 


(A, — Ar)Az, = (Az, — Ar) 4g, 


and 


(AZ — Ap)A= (AK — Ap) AK + Ap(AE — Ar). 
We have, 


IAF, — Ar) AG, III 
III(AF%, — Ar) rll 


II(Az, — Ar)IIL Pall? Az, 
(I|Prll +1) ||Arll lz, — Azll. 


IN IA 


These inequalities prove that AG “y Ar. 


Theorem 4.4.4. Let u be the solution of (4.1) and let un be the Galerkin approximation. Then, 


lim ||u — Un|| #4 [a,b] = 0. (4.24) 
n—+0 


Proof. For n enough large, 


U —U, = (Alz — AG)! las (P,U —U) + (AR - Av, 


0 —Ualle < Cy eal \IP.U —Ullg, 


4° (ak =Aypollg, + ee = Fila) 
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Numerical Examples 


Then, by theorem 4.4.1 


lima |lu— tn|larafa.o] = lima ||U — Unll x, = 0. 


Galerkin System 


By the definition of basis {e;}"_, and the duality product (4.3), the elements of system (4.14) are repre- 
sented as 


b b 
Ad) Jf Kile.tes(@ei(a) at de, iL2tjeu ete: 


b 5 OK 
Aar(i,d) = ff F*(w.thes(tei(w) dtde, 1<i,j <n, r=1,2, 


4.5 Numerical Examples 


4.5.1 Numerical tests in W*"(0, 1] 


To treat our method in this Soblev space, we give this numerical example 


Test 01 


We consider the following integro-differential equation 
1 1 
Au(ax) = | (Ie ~ tog — _ je~ t]) ul ar+ [ / |x — tlu’(t) dt + f(x), (4.25) 
0 0 


where, \ = 9 and 


ie can (252? + 38¢ + 45 —log(1 — x)(36 + 24a 4 121°)) 


4 3 8 5 4 
Es 1— a)? +92? — =a? — _24(12log(x) — 25). 
TA x+3)(l—2)? +92 a ida” \ og(x) — 25) 


The exact solution u(x) = x7, in the next table, we present the error between the exact and approximate 
solution, with the error 


erry = ||U — Un||w2.2[0,1] 
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n Galekin Kantorovich 
10 0.0750 0.0041 
50 0.0150 8.747e-04 
100 ~=0.0075 4.4641e-04 


Table 4.1: Numerical results for equation (4.25) 


Figure 4.1: Errors between exact and approximate solution of (4.25) 


0.15 1 * 
—6— Galerkin Error 
—— Kantorovich Error 
eat 
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ae 
pets 
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a 
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0.05 + 
c 
at 
\ ee eee eae va & 
0 1 i ERS tanec 2 oer eeeee Ti 1 
01 O02 03 04 O05 O06 O7 O08 o9 1 
Test 02 
We consider the following integro-differential equation 
1 1 . 
Suite) = fi Je —eu(t) dt +f le — et |Su/(t) dt + f(2), (4.26) 
0 0 
where, A = 11 and 
2 3 2 3 
f(x) = 7eMx 3(e e*)3 re e*)3 ge 1)? iG 1)3 
The exact solution u(a) = e* 
solution. 


, in the next table, we present the error between exact and approximate 


n Galekin Kantorovich 
10 0.0864 0.0219 

50 0.0172 0.0047 
100 0.0086 0.0038 


Table 4.2: Numerical results for equation (4.26) 
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4.5. Numerical Examples 


Figure 4.2: Errors between exact and approximate solution of (4.26) 
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4.5.2 Numerical tests in H'{a, b] 


To observe the error behaviour in the Sobolev space H'[1, 2], we give the following numerical example 


Test 01 


We consider the following integro-differential equation 


Vx € [1,2], Au(x) -[ u(t) a+ f “O, dt + f(z), (4.27) 


ev + et 1l+a+ et 


where, A = 2, u(x) = e* and 


i) 


e 


el 4 62 el arctan (5) 

x) =l1o + 2e” + arctan . 
a i) (a) aet 
The next table shows the numerical error between the our approximation solution and the exact solution 
which we define it as 


err, = ||u — Up||a[1,2)- 


n Galekin Kantorovich 
2 0.0287 1.9624e-04 
5  3.7704e-06  1.6025e-07 
10 6.6715e-14 ~—- 1.8394e-13 
15 1.7161e-22 9.7426e-25 


Table 4.3: Numerical results for equation (4.27) 
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Figure 4.3: Errors between exact and approximate solution of (4.27) 
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Test 02 


We consider the following integro-differential equation 


Vax € [0,1], \u (4.28) 


_ fu) yu) . 
=f yesrra’+ | yepagr ttle. 


where, A = 4, u(x) = x? and 


arctan ( 


flo) =4e? ret) 2 (vers - va¥2). 


Vat 


In the following table, we present the error between the exact and approximate solution in the sense norm 
of H*(0, 1). 


n Galekin Kantorovich 
2 2.75151e-37 3.60996e-7 


5 2.42575e-36 —5.35642e-12 
10 5.579869e-35 6.79e-20 
15 = 33.8871e-12 5.6832¢e-22 
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Figure 4.4: Errors between exact and approximate solution of (4.28) 
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Conclusion 


Integro-differential equations play an important role in many scientific areas. Through this thesis, we 
have tried to pay attention to study one specific type, which is the linear Fredholm integro-differential 
equation. We have previously explained its great importance. We tried to study the different forms that 
the equation can take. 


At first, we were interested on the equation defined in the Banach space C1[a, b]. On this basis, we 
constructed a sufficient condition to ensure the existence and the uniqueness of the solution in the Banach 
Cla, bj. After that, we built three different numerical methods: Nystrém, Collocation and Kantorovich. 
Our aim, is to compare them in terms of behaviour and convergence analysis. In addition, to search for 
the best error that converges to zero. We have provided several numerical examples which show that the 
Kantorovich method is the best of the three methods. 


In the second part of this thesis, we assumed that the kernels K; for 7 = 1,2 are weakly singular. From 
this, we were able to construct a sufficient condition that shows the existence and the uniqueness of the 
solution, which made the numerical and analytical studies coherent. Because the product integration 
method is the most common and widely used method for solving this type of equations, we applied it 
and developed a theory that explains its convergence in C1[a,b]. We did not stop there, we also built 
another method based on B-Spline functions, through which we obtained a system that contains a block 
instead of a set of blocks as in the previous methods. We also studied the convergence of the new 
numerical solution in the Banach space C1{a, b]. The advantage that distinguishes the B-spline method 
from Product integration is that we can gain in accuracy and shorter time, as the examples show that 
this method consumes half the time that the first one does. 


In the last part, we changed the space from C1[a, b] to the larger spaces W1?[a, b],p € [1,+00]. We 
were able to show the existence and uniqueness of the solution in W1 [a,b]. Then, we constructed two 
methods based on the Kantrorovich and Galerkin methods. We have studied the convergence behaviour 
of both methods through some steps, the most important of which is to show the v-convergent of the 
Kantorovich and Gakerkin operators. Our objective in this section is to answer on the following question: 
Is the Kantorovich method still the best one in changing the space? The answer is well illustrated in the 
numerical examples. 


As perspective, we would try to apply other projection methods like the Kulkurni method [39], convo- 
lution and Fourier series .... Although with this method, it is necessary to solve a system of the same 
size as Galerkin’s method, the solution obtained converges faster than the projection methods proposed 
in this thesis. We will not only focus on the linear Fredholm integro-differential equation of the form (6). 
Also, we are going to study the linear and non linear Volterra integro-differential equations. 
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